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Magnetic Polarizability of a Short Right Circular 
Conducting Cylinder’ 
T. T. Taylor 


August 1. 1960 


The magnetic polarizability tensor of a short right circular conducting cylinder is caleu- 
lated in the principal axes system with a uniform quasi-static but nonpenetrating applied 
field. One of the two distinet tensor components is derived from results already obtained in 
connection with the electric polarizability of short conducting evlinders The other is calcu- 


lated to an accuracy of four to five significant figure s for evlinders with radius to half length 


ratios ot i, . l 2 and t The sf results, when combi ( d with the corresponding results for 


the electric pol irizabilityv, are applicable to the proviem ol calculating scattering from cylin- 
ders al d to the cde sigh of artiti inl dispe rsive me lia 


1. Introduction 


This article considers the problem of finding the magnetic polarizability tensor 8,, of a 
short, that is, noninfinite in length, right cireular conducting evlinder under the assumption 
of negligible field penetration The latter condition can be realized easily with available con- 
ducting materials and fields which, although time varying, have wavelengths long compared 
with evlinder dimensions and can therefore be regarded as quasi-static. The approach of the 
present article is parallel to that emploved earlier [1] *° by the author for finding the electric 
polarizability of short conducting evlinders 

The magnetic polarizability tensor relates the induced magnetic dipole moment m, to the 


inducing field B according to the equation 


with the evlinder situated in free space. The centered coordinate system of principal axes 
used in {I in which the axis coincides with the axis of rotational SVinmetry, Is also emploved 
here In this svstem, the equal transverse components 6 and 3,, are denoted 8,,, and the longi- 
tudinal component 8.. is denoted 6 In agreement with other conventions of [1], a represents 
the radius of the evlinder; 6, the half-length; and e¢, the half-diagonal or distance from center to 
edge Rationalized NI KS. units are used throug hout 


The method described has been used successfully for calculating 8,, for evlinders with a/b 


ratios of ‘4, 's, 1,2, and 4. For the calculation of 3,,, the following relation, proved in appendix 
A, Is invoked 
la 
3,.=— 2 
A ? «€, _ 
where a,, is the transverse electric polarizability treated in [1]. It is interesting that eq. (2 


Is valid hot only lol erreulatr evlinders but lol all conducting solids ol revolution. Both 
and 8,, are newative Indicating that an artificial medium composed of short evlinder elements 
arranged ma lattice will he diamagnetic 


In the transverse case, wl ere ho simple relationship such as (2) 1s know nto exist. the prob- 
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lem which presents itself is that of finding, in the external space, a static magnetic field which 
has no normal component at the surface of the evlinder and which reduces to a uniform field at 
infinity. With such a field there is associated physically a surface current density j on the ey lin- 
der and a field free condition in the interior. This surface current is related to the field values 
at the surface by 


o] e B, 2 


where @, is the outward normal, and it gives rise to the induced macneti dipole moment 
if ) 
m= | r& j { 
As with the corresponding problem involving electric polarizability, there is no known coordi 


nate system 1n which a short evlindet can be treated by the conventional ethod ol separating 
Laplace’s equation and a different method, involving arbitrarily good approximations to thi 
current distribution function, must be emploved 

In a manner analogous to that of [1], let the evlinder be regarded merely as a geometrical] 
construct to which the surface current distribution is firmly affixed. If the uniform (applied 


component of the macnetic field Is subtracted the renaming field will be 1 ut due to the cur 
rent system on the evlinder ac ting alone: it will be characterized by a d pole hike external fis ld 
and a perfectly uniform internal field whiel Is equal and opposite to the applied field The prob 


lem Day therefore be restated In terms ol finding on a finite evi drical s Iriace with closed 


plane ends, a surface current distribution which, when acting alone, generates a uniform mag- 
netic field inh the interio! The technique employed here makes use OL this viewpoll tand cor 

sists In setting up, on the surface, a current density function which completely determined 
by a finite number of parameters, then solving for the values of these para eters so as to obtain 


maximal uniformity of the current field or, alternatively, maximal cancellation of the applied 
field, within 

The research reported here, like that related to the electric polarizab itv. Was performed 
under the direction of Prof. W. R. Smythe. The author is indebted to him for many valuabk 
discussions and to the personne] of the California Institute Ol Teeln ology Computing Center 


for their cooperation, 
2. The Transverse Problem 


In this problem, let a unilorm magcnetic field ol magnitude B by ipplied in the Positive 
y-direction. The induced surface current density j may be regarded as equal to JE where € is 
the thickness of an infinitesmal surface laver in which the current is confined and J is the 
value of the conventional current density vector averaged over the thickness of this lave 
In the quasi-static case, J is solenoidal and can be regarded as the curl of an equivalent iver- 
aged magnetization M which is also confined to the infinitesmal surface laver and which is 


normal to the surface at all points. On the side of the evlinder, let 


uy. ,£ 


. ak b)sin @, 5 
bh 
and on the upper end, 
uy. 3 0 
“t1\ p/@) Sill @» a 
B 
where F' and G are dimensionless functions as vet to be determined The configuration just 


described constitutes a nonuniform magnetic shell. that is a magnetic shell [2 p. 261] associated 


with a distributed current system 
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Evidently the surface « urrents on the side are vlven by 


ik s h Sill oO, Ss 


where the prime denotes differentiation with respect to the dimensionless argument h, The 


suriace currents on the uppel end are 


(’ pop @) S111 Or 10) 


The current distribution on thre lowe! end is equal and Opposite Lo that on the upper end. 


Typical forms for the functions involved in eqs 7) through 10) are illustrated in figure l:a 


‘) 


sketch of the ac tual line Sot current flow Is presented in figure 








a/p p/oa 
4 4 
me : [ me 
nan 
} ‘oes 
| | | + _ 
+ ——— ~ ——t—> > /b >——————— —— - 
4 
4 /| J 
/ 
A 
| _ 
+~« Ae + —— Z 
| RI Typical for for the functions describing the equiv- FraureE 2. Sketch of the current systen 
i] ; ; / ; 


The asy mptotic behav ior ol current and field heal the edge of the evlinder is determined 
bv the fact that, locally, the edge resembles a configuration consisting of two infinite perfectly 
conducting planes intersecting at right angles with the field restricted to the outer (larger) of 


the two angular domains. Such a configuration is a familiar one in two-dimensional potential 


theory In the case at hand, two mutually perpendicular components of magnetic field are 
distinguished: a) the component which is every where parallel to the edge, and (6) the com- 
ponent which bends around the edge. Each of these field components has associated with it 
a Surface current obeying eq (5). The component of magnetic field parallel to the edge is 


locally uniform in the outer space and is of no further interest in the present discussion The 
I 
surface value of the component which bends around the edge, however, increases as the edge 


is approached and is in fact proportional to / where / is the perpendicular distance from 
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the edge to the surface point in question This dependence, which is the same as that of the 
electrostatic charge densit, On a similar freely charged configuratio is obtained by means of 
the Sehwarz-Christoffel transformation The associated component of surtace current is, of 
course, parallel to the edge and evervwhere proportional to the field On the evlinder. then. 
the ¢-components of current, and hence the functions /’ and G’ must be asymptotic to con- 
stants times / as / tends to zero This. in turn, Imiposes certain requirements upon / and 


G. Quantitatively, 


I: h ~f (. h is bh: 1] 
Gipiay~f ( ,(a p is pod 12 
ad e 4 
I: ‘ h’ ~ “at h iis h: 13 
b 3 
) 
G’ pda ~ a 5 re) iis p “1 | } 
It is seen that two edge conditions, ¢ ( and ¢ ‘ €., must be satisfied The former of these 


msures CONUNUITY ot current flow ucross the edge of the eviindet the latter that the o colmM- 
ponents on side and end match one another as both tend to infinity at the edge 

The functions F and G are represented DY finite weighted SuUllS nvolving the unknown 
coefficients f,, fF». Yo; and Gms ana the y functions defined in [1] Note that the parity para- 


meters are chosen so us to make F even and G@ odd and that the parameter Whi h controls the 


singularity is given the value 2/3. 
Thus, 
\ in | ) ) 
/ * / y ( () 2") ); h 
dos Sauer 
F(2/b - = 15 
| 
and 
Jy T 2 ¥ ( = ) psa | 
( ola ad Ya 14 
t} ; 


The constant terms 


ure called the ‘‘ basi terms as in [1 Tow the they constitute a very s mpl current system 


in which the flow tukes plane ein plane rer tune@ulat loops around the eviindes the plane ot each 


loop bemg parallel to the a plane ot figure 2 Superposed Upol this Dus flow pattern is 


another pattern described by the nonbasi summed) terms Ld ina lo The addition 
of this pattern gives a gentle curvature to the current lines and introduces the necessary singulat 
behaviol nt the edges ( learly the first eda condition demands Lhiat j and henceforth 
the single symbol /, will be used as the coefficient of the basic term in bot! Fand G 

In relating the self field of the current system to Mé nh macnetl Cruial potential defined 
such that B VU, will be used Th 


duals of the corresponding boundary conditions for the electri pote ntint ata chare@ed Suriace 


e boundary conditions at the magnetic shell are the 


that Is, the normal grad ve nts on the two sides of the sii || nie equal Dut the potential (LUES 


differ by the local intensity of the magnetization. Quantitatively, if | and Il are two adjoining 
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regions separated by a magnetic shell, and if @, is a unit vector pointing from I into II, then 


',, ol 
() 19 

On OQ? 
Uiu—Ui=—p.M,¢ 2) 


Application of these boundary conditions ana of other well known procedures of potential 
theory leads to the following expression for the pot ntial in the interior of the evlinder due to 


the current S\ stem on the side 


.D, ( 2 , , , r 
E : | F(z b) cos f d a; kaj] kp cos kz dk sin 47) 21 
ch ( a - os : 





The Fourie! cosine transform ol the I function, enclosed in the square brackets above, Is easily 
expressed in terms of Bessel functions with the aid of eq (A9) in [1]. The function /;(kp) cos 
} Sin oO, which is revular In a neighborhood of the origin, may be represented nm a spherical 


harmoni expansion nbout the origin Thus 21) becomes 


/ r 60 ab y ; r P g kb \ I | "J » kb) 
cB ee | i ebyv2 + 2 J kb ] 


a 


hak ka) (ke)*? A} 2 (= ¥ P3 cos #) sin o. 22) 


Here and elsewhere, the singularity parameter v’ for the nonbasic terms Is understood to be 
equal to 2/3. The spherical harmonic expansion of (22) converges for all r<a 
With reference to the ends of the evlinder, the interior potential venerated by the current 


distributions on the two ends ts given by 


/ pP.Q, 2) a 


. l ¢ : ; 
B | G p’ a)eJ kp’ p'dp" | ka aa” kp) cosh k dk sin @. 233) 
CD Cc a 


This, in turn, may also be expressed as an expansion in spherical harmonics which converges 
for all» h 


/ r 4.0 ~~ a [ i. J ha \ ef “1 . | m ka) 
cb am a ha => ka 


] <1 ig 
hae~*’ (ke dt aN ( - ) < cos 7) sin @. (24) 


Evidently the sum of (22) and (24) constitutes an expansion of the interior potential due 
to the total current distribution, convergent within the largest sphere which can be inscribed 


in the evlinder Such an expansion may also be expressed in the form 
\ \ r\* ; 
» > ZZ - S* {2 t 2° + y > d Zs |( ) cos 6) sin o, 25 
{ 


where the Z’s are obtained from (22) and (24) by integrating [3, p. 137] over k. These Z’s 
become the matrix elements in a system of simultaneous equations which may be solved for 
the f,. f,.. and g,, such that the coefficient of (r/¢ cos @) sin @ will be unity for p=0 and 
zero for as many p70 as possible The second edge condition, also included in the system ol 


simultaneous equations, 1s simply 


“ 
a“ 


. a . . 
> Sat{ j ) > g,,=0. 26) 
; 








If one makes use of the fact that kak, (ka kak, (ka A (ka), the expression for the 
matrix element Zs becomes 


ab/2 ; J eh 
ZL ( - ) ap 5 ; kak ha A hea Lh + /] oy 


i . 


Integration vields the following result 


The matrix element 7;" is simplet 


( 2p 2). ha 
This becomes 
Zi" 1)"2° Tm p 5 (Cm D-+- 2 °(“) 
r(3)P (2p+-3)P (2m+3-49’ oN 


The basic miatrix elements Z and Z are obtained by setting DOT ind . equal to 
zero in the preceding formulas 
Thus, 


, 1) ?2° CP (p+ 3) 0 (p+ 4) a*b 
a r(3) 0 (2p+-3)T ( | (y ,) 


y 2 


/ 4 
[ p | 

The second term in (33) vanishes unless iL 0, in which case the nypergeometric tunetlon is 
easily summed 

| I h h f 

) 
re ade ty 42) ’ 

~ a= ( ad 

Let the svimbol rs) te equal to uUnItY lor Pp () and to zero tor p “~(): ther 


In (25), the total coefficient of f, is the sum of Z, and LZ, . This sum ts denoted by x and is 
given by 


; 1)?2 (p+) (p+) a’b | ; 3 6 


2p+3 2 2’ « 
hi (v3; 32)" (nth (vps 2) 
Using the third Gauss relationship, as listed in [4, p. 103]: 
ys 1)?2”T'(p rp 4) ah , (oak a 1-3. 0) (ag 
r(s)r(2p+3)re ( Po a’ € 


This has an even simpler iorm, obtained by an elementary transformation ol the hypergeo- 


metric function and a simplification of the vamma tunction multiplier : 


: 1)?r(p+4) b b: 
ZL 6 var / = I ( p ty» —P; 33 ): (39 
2r(4) P(p+2) ¢ : 2 
The three formulas, (28), (30), and (39), contain the information necessary for a numerical 


calculation of the matrix elements. If both the electric and the magnetic polarizabilities are 
being calculated lor the same evlinder, much labor can be sitved by comparing the expressions 
for the matrix elements in this paper with those in [1] and exploiting the similarities. 


The magnetic dipole moment, given by eq (4), has only a y-component, 


: ; . 
Mi - ] SIN @—d) COS @O adod + bi Joe COS D Joe SAD o)ddpdp 1() 


Using expressions (7) through (10 


uy mm vy ’ a , “(ar = ” 
5 | ( I Sin” @ ak cos* @ ) adgd: + b ( ——(1 cos” v8) G sin® ¢ ) dépdp; (41) 


. . JO. , p 


- 7 Trad } | |- (uk I du Gy uG’ du (42) 
Ly mM , ee 
BR rab luk 2| kdu uG }. (43 


The first edge condition requires that Fd Gi) hence 


By parts 


Uy Ih 


B Pnra*h du (44 


The quantity hy B is the transverse polarizability and the ceometrical volume M is easily 
recognized The orthogonality relationships frivenh in 1] aid in the integration of F. vielding 


the following formula for the polarizability 7 


(1 
unBe: | fo ri | (45) 
Dp , ehe eC 
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3. Results 


The LWwo distinet tensor components of the magnets polarizability are given mn table - 


Actual values of the expansion coefficients for the F and G functions, which describe the 


surface current distribution via the equivalent magnetic shell for the transverse case, are 


included inh appendix B. These data, upon which the 6 results of table 1 are based were 


obtained by solving sets of 18 simultaneous equations Experimentation with fewer equations 
has led the author to the conclusion that the accuracy is not quite as good as that obtained 
for the electric polarizability in [1] In other words, the 8,, figures have a probable error of 
about 2 units in the least significant digit for a/b equal to unity and about 10 units in this 


digit for a h equal to 4 or to ‘ The OD fivures, on the other hand were obtained irom the 


electric polarizabilities simply by the application ol eq 2?) and theretore ehyo\ the Same 
accuracy belonging to the latter 
PABLI 
p ; 
‘ . 


As an additional check upon the numerical calculations in the transverse problem, a 


procedure analogous to that emploved ith L| Was used Consider the total field B which 1s 


the sum of the applied field and the field due to the current S\ stem on the cy lindet \ normal 


component B,, of the total field at any point on the surface of the cylinder indicates that an 


error 1s present If represents distance uway from the surface in the outward normal direction 
at a particular surface point, one may calculate the displacement of the surface An which, con- 


sidering the local value of the normal field gradient, would reduce the normal field component 


to zero In other words, let the following equation be solved for A 


an n=O 16) 
' On 
This An may then be compared with the gross dimensions of the cylinder. The local field 


gradient 0B,,/0n depends upon the current distribution at the point question and may be 


calculated by appealing to eq | and Lo thre solenoidal charac ter of B To first orde} 


has for any point on the side of the evlinde 
Ob lo 0) ) 17 
op ado: o ' 
ob B sin #8) ad 


>. = (/ 5 FF? ) 18 


Let the p-commponent of the field due to the ¢ urrent svstem alone be denoted B Then 
B, B B Sin O, A? Ap, and one obtains tor the latte) 


il i 


b 


I 
Ap BB sin @ 
1 
ad ad 
F'+;, I 
h 


The quantities 4,, Fy F’” and finally Ap/a were calculated for 0, that is for points on the 
“equatol of the evlinder: hence the term “equator check”’ In appendix B. 


The same method applied to the upper end of the cylinder vields 


ob ise jee , 109 
wl 11 Qin 


ov) 
O Op 0 Oo 
Ob, Bsing o .. « 
( G ): 51 
O ad / 7 
where &@ represents p/d Substitution of the expansion 4,44 Aew Au ... for G. which 
s known to be an odd funetion cIVeS 
ob B Slli @ 
S.1,u— 24.150 52 
O ad 
Since both sides of (52) vanish at p=0, that ts at the pole ’ of the evlinder, one must differ 
inte wit respect to pih order to obtain somethin oO useful Thus. 
J Bb sin 4) i/3 sili @ —— ; 
Sal G () 52 
OpO a ¥ 


‘urthermMore iol wri pont ol the lppel end 2B - simply B the field due to the current 
system alone Finally 


A od B 


O , 
h 1ha’’’ (0 >. | B sin rt 54 


Values of Oo Op. G’”’’, and A bh were cale ulated lor the pole; data are given under the heading 
ol pol cheek” in) appendix B Generally speaking, the pole checks are satisfactory but are 


not as rood us the corresponding chee ks in 1] 


4. Appendix A: Relation Between the Longitudinal Magnetic and the Trans- 
verse Electric Polarizabilities of a Conducting Solid of Revolution 


} 
> 


Let figure 3 illustrate a tvpieal solid of revolution formed bv rotating the curve p= o(2 


itootit thre AXIS 
+ 
e. 
> 
Fs 
\ a 
+x 
| R} ( . pica Jed a 
} ¥ j 4 si | 
Y 
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Let the unit vectors @, and e, point along the tangent and the outward normal of this curve, 
The variable s, associated with @,, measures distance along the curve and encom- 
S. It is sufficient 1 has a uniquely 


respectively 
passes the solid as it varies over the interval 0<s 
defined tangent and is normal at almost every point in this interval 
If the applied fields consist of a uniform magnetic field B in the z-direction and a uniform 
electric field / in the x-direction, the following is to be proved 
My Mt l Pr 
> . ’ \] 
b Z2eok 
In the electric case, the scalar potential satisfies the conditions 
V Hey p, 2) cos @; . 
f,=0 on the surface of the solid 
> \2) 
f;~p at infinity; 





4 


v°7V—0. 
The magnetic case, on the other hand, is analyzed with the aid of the vector potential which, 
This component must vanish at s=0, s=S, and at all 


by symmetry, has only a ¢-component. 
intermediate surface points; otherwise there would be a normal component ol magnetic field 


at the surface. One finds that 
B ) 
A—e, - f p 
f.=0 on the surface of the solid; ‘ AQ 
f,~p at infinity; 





VA (), J 


Application of the Laplacian operator to either cos @ 7; or to @g fo \ ieclds the same equation 1n p 


In other words, 
(A4) 


oO 1 Oo O° / 0 
e” * ef Je 


and 2. 


Op” p Op 
Since /f; and f, obey the same differential equation, have the same boundary conditions at the 
f= 


surface of the solid, and have identical asymptotic behavior at infinity, it is clear that / 


The surface charge density in the electric case is 
o 6 @, -VV=eF cos ® (@,-VJ) 15) 

and it is easily shown that the z-component of the induced electric dipole moment 1s 

Dr=€o fix (e, - Tf) p*ds. Ab) 
The surface current in the magnetic case as given by (3) is 
» 9 B . 
J e, X(V XA) =—— €, X(V XCof) AZ) 
Lu ~h 


0), this reduces to the following 


On the surface, where 
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and the magnetic dipole moment is obviously 


Br 
m,=- e”- Uf) p*ds (A9) 


A comparison of (A9) with (A6) shows that the asserted relation (A1) is true 


5. Appendix B: Values of the Expansion Coefficients for the Functions F and G 


Values of the expansion coefficients for the functions F and G and the results of the checking 
procedure at pole and equator are given in tables 2 through 6. The notation “A( p)” means 
‘A times 10?.”’ 








TARLE 2 TARLE 4 
l 
+0. 13640434 (4 basic +-(), 12092027 
+0). Y4335081 +) 0. 37183474 (+0 +-0. 54323936 (+0 0, 29402467 
+0. 47094227 (+0 IS5H3102 (+0 +-(). 92008198 l 0. 12410492 
+0. 15962976 (+-6 0. 97767116 l +0. 17619634 ] 0. SORGS2RS 
+() l +-(), 38217146 2 0. SYQ453649 
l +0. 67793221 3 0. 44063915 
55233089 2 0. 6676957 4 0. 29238560 
+0. 16203904 (—2 0. 16182197 
+0. 41694492 3 0. 68464716 
+0. GOL3R919 4 0. 19415309 
+0. 15575188 (—4 0. 27441144 
+0). DOO4S427 
+0). 17022684 ( 
1319103 (—8 /Bsiné 0. 99999908 
Ap/a : neglibl 
f B n ee FHsin 4 —0. OOOOO4SS 
\é 
Azit 0. OOOOLSI2 
R,/f I +0. 11001902 
\ +). OLV464905 
TABLE 5 
TABLE 3 
). 12407197 l 
+0). 33590275 0. 25652350 
+). 30150397 ] 0. 11210208 
13390263 +0). 28485520 2 —). 72413481 
0. 523863908 
0 238 (+0 0. 32481666 (+0 0. 37934754 
) 2 +) 0. 14321008 +-() 
0 254 l (0). 97128304 l 0. 25606044 
). 1676487 l 0). 67455522 l 0. 15530579 
+0. SOOG2 106 2 0. 36913235 l 0. SO562157 
0. 344986033 
+0. 13805064 9 0. 11074155 | 0. 11632626 
+0). 32120205 ry 
6756068 { 0. 28842547 
0 ;WHTHOOL 0. 46632080 
{ 17047 ( 0. 36819245 
Ry sin @ 1 QOOOO2SS 
R,/RBsin ¢ 1. QOOO0004 
\ me t \p/a +0 00000 L15 
/ / 0. OO68 1408 R, sin ¢ 0.000000 LS 
\ (2256388 Aci 0.00000 LOO 
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Accuracy of Monte Carlo Methods in Computing Finite 
Markov Chains 


N. W. Bazley and P. J. Davis 


\Iav 20, 1960 


boxpe ments are icle Vil tthe Marko I presented by the children’s game of 
( es nnd | idade rs Statistics suctl im tle iveruge engt! oft play nre computed on the 
IBA 704 fro 2 1 " ted plavs of the game These Monte Carlo results are then com 
he “exact itio obtained by powe the matrix of transition probabilities 
( ~ t i | \ 
l. Introduction Since the plaving piece can never remain at either 


the top of a chute or the bottom of a ladder, there 
since Monte Carlo methods are generally applied re only SI squares at which it can stop. 


tliat are too difficult lol othe types ol The plaving prece Is off the board nt the Start of 


0 problems 
there is seldom OPporvunity to colpare the game and the plaver tosses the die so as to deter- 


nalvsis 
to 6 He then moves his 


elr results With “exnet solutions In this papel mine a number from ] 
yperin eChts are described with thre Markov cham play Wie pre ec through the Ssanie numbet ol squares as 
by the ehildren’s game of Chutes and the number on the die. Should his playing piece stop 


resented 


adders I particular, the ecaleulation of the aver 7 | the bottom of a ladder he must unmediately move 
re lengt! a Ol a play ol this Yale is considered to the top ol the ladder He then repeats the throw 
in exact value for [is found according to Markov of the die and moves according] If at any time 

aE : : — he lands at the top ol a chute he must proceed directly 
hain theory Phe \lonte ( arlo eal Wation ior i to its bottom wr ae thee nent din hoon The vame 4 
volves 2 stnulated plavs of the Cule Further- comeledad when the plaving piece reaches square 100. 
ore, thie probabilities ol being In a given state by The hoard actually has two winning positions, 
fixed number of throws are computed theoretically since there is a ladder from square 80 to square 100. 
OG CAPCs tt ntally One requirement Is that square LOO he reached by an 


It is shown that well-known theorems of statistics 
ply that the Monte Carlo results should converge 
o the true solutions by ‘oA i laws Numerical 
ilues are shown to follow these laws All work was 3. Theoretical Solution 
ried out on the IBM. 704 compute! 


exact throw of the die. 


The game can be interpreted as a finite Markoy 


2. Description of the Game chain with 82 states. These states are indexed by 

o, 3. . Sl, where the state 7=0 corresponds 

First the game of Chutes and Ladders is deseribed to the starting position and the state 781 corre- 
detail The game is considered to be plaved in its sponds to the plavi ig piece being at square LOO All 
solitaire form The game involves a plaving piece, squares except those at the top of a chute or the 


me die, and a playing board which has squares | bottom of a ladder have a corresponding state in 
umbered from 1 to 100. Certain pairs of these 100 | the chain 


squares are connected by ‘chutes’? and others are Let .Y,(n), s==f, 1, . 1» Sl represent the prob- 
onnected by “ladders ia The purpose ol the chutes ability ol being in state i after i throws of the die. 
sto delav the game by foremge the plaver to return Then  Y,(0 i and N,(n)=—0 for n>0 Also, 
to the square al the bottom of a chute if he should ‘ b : : 

: Ae 1 for every 220, Let Pp be the condi- 


ind on the square at its top On the other hand, the hat 
ulders speed up the game by permitting a plaver to 


the bottom of a | “ona! probability of moving from square 7 to square 7 


divance directly from the square a , 
: on one die toss whenever square 7 1s occupied, We 


ladder to the square at its top, an 
squares in between There are LO chutes connecting 


the pairs of numbered squares (16,6), (47,26), (49,11), ; : 
56,53), (62,19), (64,60), (87,24), (93,73), (95,75), Ek rob of a 


l 
| thus bypass the 
see immediately that 


ind (9OS.7S However, nine ladders connect the J on nh Losses | 
squares Las $.14), 9,31), 21,42), 28,84) Pro] f tr it 
6.44). 51 67 71.91). and 80.100 & Prob 0 rans lon Prob of state P = 
from state 7 to . » (29) 
— ae _— on 7 l tosses 
( 1 ( State 7 on one toss 
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In other words, the elements .Y,(n satisfy the 
relations 
Ps () S| 
X,(n)=D py XAn 3.2 
} Lue 
() 1<4 S|] 
where .Y,(0)= If we let NG be a 
l, } 0 
row vector whose components are | TL Wwe ean 


>”) 


write (3.2) in the form 


X(n)=X(n—1)} 


X(0) the vector 
(0), —- * 


where 


X,(0), X; 


Is 


0) and / 


0 l 2 } } 5 6 7 

0 1/6 1/6 1/6 1/6 
l 1/6 1/6 1/6 1/6 1/6 
2 1/6 1/6 1/6 1/6 
} 1/6 1/6 1/6 1 
1 1/6 1/6 1 
5 16 1 
6 l 
7 1/6 

Ss 1/6 

9 1/6 

10 1/6 

11 i/6 

12 1/6 

13 

14 

t 

66 

67 

OS 

69 

70 

71 

79) 

74 

40 

76 

78 

79 

SU 

81 


W hose 


components are 


"is the nth power ol 


TABLE 1. Partial listing 
S ’ 10 1! 12 13 l 
1/6 
1/6 
1/6 
6 1/6 
6 1/6 1/6 
6 1/6 1/6 1/6 
6 1/6 1/6 1/6 1/6 
1/6 1/6 1/6 1/6 1/6 
1/6 1/6 1/6 1/6 1/6 
1/6 1/6 1/6 1/6 1 
1/6 1/6 1/6 1 
1/6 1/6 1 
1/6 1 
l 
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p 
pP 


Psior P 


The matrix elements p,; will be called the transitioy 
probabilities; certain of them are given in table a 


Ps 


We leneth, LZ. 


conside! 


now the average of thy 

came This may be computed from the expression 
b=)  nlXa (1 NXg (n—1)], 3.4 

where the term XY } \ # iT | represents the 


probability of reaching the final square on lie Toss 
number 7 

\n alternative expression which can 
computing ZL is given as follows. 
truncated matrix 2? 
last and the last 
matrix 


be used fo 
Kirst introduce th 
formed by crossing out both the 
row column of P? Then form the 


r 3.5 


l n 
rt 69 70 71 72 7 75 6 77 78 79 80 81 
6 
( 
' 
a 
tr 
6 1/6 1/6 
1/6 1/6 1/6 I/€ 
1/6 1/6 6 1/¢ 
t ft) ' ' 
{ Lif f 1/6 
{ 6 1/6 1/6 
6 1/6 1/6 1/6 
f 1/6 1/6 1/6 
1/6 1/6 1/6 1/6 
1/6 1/6 1/6 1/6 
1/6 1/6 1/6 

1/6 1/6 1/6 V6 

1/3 1/6 1/6 1/6 

1/2 1/6 1/6 

5/6 1/6 

l 


Sition 


e ], 


rf the 


CSSIOT 


S thy 


toss 


d fo 
e thi 
h the 
n the 


6 1/6 
6 1/6 
6 1/6 
6 1/6 


[t is shown on p. 329 of Kemeny, 
that Z is given by 


with elements 7 
Snell, and | hompson , 


L > T, , 3.6 


which is simply the sum of the first row of the matrix 

I-P; 

Equation (5.5 that a direct numerical 
solution by successive matrix vector multiplications 
can be calculate the vectors .Y(n) This 
solution may not be feasible for full matrices or for 
chains with a large number of states, since powering 
the matrix may be subject to space and time limita- 
tions and roundoff error. However, for Chutes 
and Ladders only 476 of the 6,724 transition prob- 
abilities do not vanish. The XY i were easily cal- 
culated on the IBM computer by taking the sums in 
,9) over the nonvanishing elements. The time 
required for each matrix vector multiplication was 
about 2 sec. To illustrate the results we have 
presented the vectors (7) and X(31 in table 2. 
The vector .\(7) is of interest since the computations 
showed that uN. (nm (0) that the game 
cannot be completed in fewer than seven throws; the 


shows 


used Lo 


for Vi 7. SO 


» r " P » 
vector V(31) is of interest since .Yy,(31) has also 
been computed by Monte Carlo. 
PARLE 2 The 0 \ ” \ / 
\ \ \ \ 
) 20 0 
{ SN 706. 10 
1 t 7413 10 
{ 12800 1 
Ww 
tit if 1 
{ mye | 7 41 ] 
s sie ts » sU44 lf 
. RIN { q 3408 & 11 
$518 4 { SN 2 
{ wl 4 SWS 
14s tty SUS i321 lt 
‘ ) yy $21] 1 
~ 4 » y 486) 
- tte! : 
{ j 15.500 X 1 
WIS s 7] 
‘ re OH 
j S00 ) 658 
s 1466S 
‘ s yaa | 
j } s4 y. 
; ss 5 i 
} SS 4 tO $504 i 
| s { lt 
’ ( ) s 
ss ; l 
‘ s H10y is ry s l ] 
) xT ) | If “4 lf 
i 14 $208] si4 l 
j S7S 17404 St x If 
+ i wt sx i 
“ OM { st Hx 10 
ve j SAS 87243 1 
{ “787 ] x If 
x 2s Tl s4 l SOOS6H XX 10) 
) ‘ ix li SM 7YLISXI10 
S444 7s s 1034 x1 W411 X10 
) 29304 X10 74177 i ) H4u32 x 4428 x 10 
» ( ~“ } ) 44su2 & If 
ts 1x {8004 
= Is G. Tl psor An t 
I e-l I York, N.Y., 19 


The number L was computed from eq (3.5) and 


uve 


as the average length of the game, The same result 
was obtained by summing the series in (3.6). The 
machine computing time was 12 min for the former 
calculation and 15 min for the latter. 


4. Simulations 


The random element in the game of Chutes and 
Ladders arises from the toss of a die at each move. 
This was simulated in our work by employing the 
pseudo-random number generator of Taussky and 
Todd That IS, a Sequence of pseudo-random num- 
bers ¢,, te, where 0<t,;<1, were generated by 
the rules 


are 2 ] 
(4.1) 


ry (mod 2°), r,=—1. 


The unit interval was divided ito six equal sub- 
intervals and the result of the jth die throw was iden- 
tified with the subinterval (first, second, , Sixth) 
in which ¢, fell. The game was then played N times 
calculator. The total number of throws 
necessary to play these N games was recorded as 
be, LL, where L; denotes the length of each game. 
The average length of the game was then computed 
as (1/N)>°*_,L, Furthermore, the fraction of these 


eames which finished in 31 or less throws was com- 


on the 


puted. For N=2"* the computing time was about 
15 min. The results of these calculations are viven 


in table » 


TARLI > Res ilts of the Vonte C'arlo calculations 
l af Fraction of 
\ — 7 I games 
ny To finished 
: $1 tl ‘ 
7 469 
1) VAY Tala 
s $6. 117 55469 
39. 480 454 
38. 701 i414 
i $Y l 47168 
48 44. O32 4007. 
Hurt 1, USS i7u74 
xa? 39. 469 47603 
‘s4 3y. 40 i174" 
EX \ 9. 224 $5004 


5. Convergence 


We first show how (] N) > Stel should converge to 


L=39.22. We introduce the distribution function 
f(n NXg; (n)—Xg,(n—1). 5.1) 
O. Taussky and J. Todd, Symposium on Monte Carlo n john Wiley 

& Sor New York, N.Y., 1956 








The mean ot /(7) 1s simply the uverage leneth oft the 
fume L as viven by eq 3.4 Furthermore. the 
standard deviation for the distribution function 
is given by 

S* (n—L)*[ Xe, (0 X 5.2 


The \ alue of o for ( ‘hutes and Ladders Wis ile ulate a 


( 


according to (5.2) with L—39.224 and the known 
values of \ / \ / | We found 2) 222 
Let Fa he. Lin be regarded as a random 


sample of size N drawn from the population with 


distribution (7 whe re the moment renerating Tune 
tion of f(n) exists It is well known that eve 
though /(n) is not normally distributed, the “sample 


mean” 1/N> 0° L 
which approaches a normal distribution with mea 


\ Vas .V becomes infinite 


has il population distribution 


Land standard deviation o 


Hence the probability is approximately .95 that the 
average leneth ol play computed from the first .\ 
simulations lies in an interval centered at 4 wit! 
leneth lo y.V. This bound involves which has 


heen computed in the case of Chutes and Ladders 








In figure 1 the theoretical values of / na ur 
used to plot the interval (4 2a/~x ies 2a/y.N) as 
afunction of V. The values of (1)N)>3>_,4, are also 
plotted and it is seen that they lie within our confi 
dence interval 

Figure 1. Co ve 


Now consider the probability ol completing the 
raume by a specified number of moves Si \ >| | 
Pp Rats then the first .V simulations can be rt 
carded us N Bernoulli trials for Ww hich the probability 
The probability ol 


/ 
I 


of success in a single trial is p 
exactly v successes is given by the binomial distribu 
tion 


‘ ] i} 

' \ ' \ ys \ ) 
Kor laro \ now ( ( >] \ 7 s ICCeESsag 
Is ipproximately orm distributed wi meg 
p and standard ce \/ \ \pproy 
mately Qs percent | re (if ~ normal \ 
curve lies thy rite 

Zy i . Zu) 
\ \ \ \ 
he st 95 peres ( ( 
\ 5 
\-\ \-\ 
Thus j/' eT > Ih) tl \ 5 is = (| 
figure 2 
ae 
42} 
} 
| 
RED. ¢ 
6. Existence of the Moments of f(n) 

In the previous se¢ »] Was assumed at tl 
moments >‘ } \ \ eCXISL Tol ill f 
If we expand \ terms of the characterist 
values ol the matrix lO! exampli See \lontroll 
p. 419) it follows that these moments will exist if / 
has A l as the only ( Value with \ How- 
ever, our Markov chan has one absorbing stat 
which can be reached from at other state The 
probability of reaching ¢ absorbing state is inde- 


Dut t} 


necessil 


that 


iis 
Crantmache! 


Is 


pendent of the initial st 


and sufficient condition 4.5 


1) 
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ind the convergence according to the 


there exist only one eigenvalue with magnitude | calculations 
V). Thus, the moments exist; also, the average length of \ ‘laws is evident from figures 1 and 2. It is 
» e gauime can be defined by eq ».0 worth noting that the ‘“A ? convergence implies 
~— that if an additional significant figure is desired in a 
| 7. Concluding Remarks Monte Carlo calculation the computing machine 
ane Lime requirements ure increased by it facto! of LOO 
oa The \lonte Carlo calculations are in agreement 
rma! th tty heoretical results obtained from the matrix Pape 64B4—36 
7 

ed 

th 

}- 

‘ist 

‘oll 

if P 

[ow 

tal 
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rice 

Sul 

tha 
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Error Bounds in the Rayleigh-Ritz Approximation 
of Eigenvectors 


H. F. Weinberger 


June 


The difference between 
ipproximation wu 
Lhe 


anv eigenvector 
is bounded in terms of the 
Ravleigh-Ritz upper bounds x The 


ipproaches zero with « \ 


1. Introduction 


The most common method of approximating the 
eigenvalues Ay <A, of asymmetric linear oper- 
tor A is the Rayleigh-Ritz method [1,5,19]... This 
reduces aun eigenvalue problem on vw space of i large 
or even infinite number of dimensions to an eigen- 
value problem on a space of relatively few dimensions. 

If the desired eigenvalues A, are characterized as 

inima, the Rayleigh-Ritz approximations x, give 
Ipper bounds for them 

Along with the upper bounds «,; for the eigenvalues 
the Ravleigh-Ritz method yields associated vectors 

It is to be expected that these vectors approxi- 
Further- 
ore, Ibis to be expected that the better the eigenvalue 
A; 1S approximated by Ki, the better will be the 
pproxXimation of w, Lo the corresponding eigen- 
ector. 


ate the eigenvectors of A in some sense 


Indeed, this is easily seen in the case of the first 


genvector. If the unit vector w, is expanded in 
terms of the normalized eigenvectors u, of \, we have 
uw Ja Us. | | 
i=! 1.2 
-_— 
a 
Sova K . 


Subtracting \, times (1.2) from (1.3 


MSA we find that 
kK 
ad a | } 
i Ao—A 


r equivalently 


“a SS a 
1S (wimg ms) | 1s 


In general the A, are 


Figure icket ite the literature references at the end of tl 


unknown, and we must express 


differences between the eigenvalues A 
bound for the 


. 1960 


of a linear operator A and its Rayleigh-Ritz 
of A and 
» and w, 


difference between wu 


our results in terms of the «x; and any lower bounds 
X, for the \, that may be available. Such bounds 
can be obtained by various methods (see for example 
| ,2,3,4,5,6,10,11,12,13,14,19,20,21,22,23,24 25, 26]). 


If \,>,, the inequality (1.5) leads to 


; i—X 
[ 5 (w Uy. W, Uy ] > | ~ — (1.6) 

2 Ao —A 
This inequality shows that if the maximum error 
K \, is small compared with the interval XA hi, the 


uv, is smallin norm. The bound (1.6) 
is sharp in the sense that equality is attained when 
A x, and a, 
trivial for «x, >Xz 
In this paper we generalize the bound (1.6). 
vive a bound for the norm of w, 
fIVen bounds iN and X ; 


difference w, 
0 for 7>2. The inequality (1.6) is 
We 


p—U, in terms of the 
This bound is again sharp 


in the sense that equality may be attained. The 
bound for (w,—uU,,w,—U,) Is small if the maximum 
error xk, X, 1s small relative to both r, Ky and 
— It becomes trivial if «,>X,., or x», >Xp>. 


The case p=1 gives an improved but more compli- 
cated version of (1.6) 

If X, is multiple eigenvalue, we can only expect 
w, to approximate one of the associated eigenvectors, 
Hence if \, lies near to several other eigenvalues we 
must expect w, to approximate not wu, but a linear 
combination of the eigenvectors corresponding to 


the nearby eigenvalues. This approximation is 
established in section 3. 

It Is possible to find 2 bound for (Wy, Uy WwW; Uy) 
by determining to what extent w, satisfies the 
eigenvalue equation. Such bounds, which involve 


Aw,—k»W, Aw,—K,w,) have been found by several 
authors {12, 20, 23]. Our bound, however, involves 
only the A and the lower bounds 4 It should be 
particularly useful in the case of differential opera- 
tors where the Rayleigh-Ritz trial functions may 
not be sufficiently differentiable to give a finite 
value of (Aw, Aw,). Our bounds are established 
by algvebraic means for the case when A is an Nx N 


matrix. They are independent of V. Consequent- 
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hold for 
first p 


ly, the bounds also infinite-dimensional 
operators A 
vectors are approximated uniformly by 
sequence c1, of \ This is ee rtainly 
the case if Al its completely tinuous It 
holds under the weake that 4 
discrete eigenvalues defined by a minimun rieaaNi 
mum principle TITAS 


spectrum Schroedinge) 


W hose elvenvalues rnd eiwen 
those ol on 
matrices 

Col ilso 
condition have p 
These nitust lie below 


Thus uA 


operatol corresponding 


Cont 
WOUS aitia’ I< il 
to both bound and unbound 
States 


The 


values of thi 


Tact that thre wenvalues \ nie Stattonary 
Ray leigl 
approximation ot the 


quotient l¢ nas tO thnk 


eigenvectors Worse than that 


of the eigenvalues In fact, the bound (2.4 shows 
that the square ot thre norm of the error ha j S 
ol thre order » \ 

The erro! bounds In sections 2 and 3 are im the 


sense of the norm If Ais a differential operator 
functions It is 
tO approximate tle 
| particulal pot An adaptatior ol the method 


of Diaz and CGreenbere 7. Of which le; 


its eigenvectors wu ure often of 


interest value of the function at 


DOLIEWISE ition is presel ted 


uppronril 


2. Separated Eigenvalues 


Let cf be an hermitian JN \ trix It is 
lineal oOperatol On Euclidean \ Space Li I eat 


usual scalar product between 
on this space be denoted by 7) 

Let the eigenvalues of A be X N \ 
The corresponding eige) 


7 My 


vectors nie adel 


l 


so that the « quations 


Au NU 


ure sul isfied 
Then (2.1 niplies 
Aus u Ad \ +S 


Rit 


bo Mas 


nethod 


Ravi Bae 


We suppose that thie 
$ Is applied to find uppel 


} 


I, 9, 19} 


a .™ o # \/ 4 


for the first \/ V eigenvalues of 4 
by choosing \/ linearly independent 
Ms ry and finding the roots 

of the secular equation 


This Is done 
vectors 


det Ave Us A\fla,t8 } a@.O | 7 \/ 4 5 


(Associated with each of the «, is a linear combination 


ot unit le neth such that thre 


\/ lineal equations 
Aw , / ) | \/ ) 


} 
| ( 


are satished 


ins that w i ist Sore (“en of 
degree ol approximation of the to the eigenvalues 
We seek fo deter e¢ iro (his degre Ol approx 
mation of the x, to the venvalues the degree of a 
proximation ofa part ir Ravleigh-Ritz eigenve 
mw, to the eivgenvecto of | 
We assume that A \ A and that 
bounds \ ana A il sufficiently rood 
\ ) 
and 
\ | 


Quy problem: is thi following 


\/ vector! 


Criven the 


wm My, Satisivu 2.0 find the largest valu 
of the deviation (9 / of w, from the p 
normalized eigenvector of anv matrix of satisfvir 
the eq 2.9) and having ¢ ve) values \, satisfvi 

2.10 


If the eigenvalues \, of A are given, 1 Is complete 


normalized eigenvectors 


specified hy prescribi : _ 
These 
thew, ce 


Let 


eigenvectors ‘,. form a basis, 


so thi 


! he 


written as 


near combinations of ther 


sSimce thre /,; Sactist\ Soi Lilie 2.0 the eqs (2.8 
and 2 Y ecome 
* ’ » 14 
ana 
2; Na ad ) | \V/ 2.1 
On the other hand 
uf / / 2 ] ad it 
If uw, satisfies (2.2) a 2.3), so does l i} 
choose thes rl) to make a non-negative Thi Ss, Uv 
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problem of mi 
to that ol 


xiImizing (uv u 7h u is reduce 
minimizing inder the 


, and (2.10 


constramts 


2.14), (2.15 


Hs 
We first keep the eigenvalues X, fixed and unequal A r ~ ‘(a | mS X, hk ‘a 
ys! [fai is the minimizing set of coefficients, we find by Bacon 5; , 8 


lirect differentiation that the equations 


VJ 0. / =p 2 IH 
LS* (r. d.8.e)d (0) Po \ Ba 
2% a 
cf eee \/ Suppose how that exactly L of the coefficients a2 
are nonzero 
2g ist bye satisfied The ros and Nos “are Lagrange mul- 
pliers, and are svmmetri a3, =v i) B=P,, Bo. ...5f 
Ol] > 5 » 7 
4] , 9 18 re) ih 3 wae 
By 2.24) one of the Bg Pp The term in braces in 
9 S ‘Ba 2.26) times the product of the (A K3 with 6,#p 
4h | ‘ - Isa p lv al eree Lm A,. Hence if vanishes 
lo solve the eq (2.17), we multiply by @) and sum Ap deg of degree / 
te th respect. to [sino (2.14) and (2.18) we find for at most 1 values of 7. Consequently, a) #0 for 
ttl CS *( Si1Y ys iti ys 2) c ( — . 
ues ) at L<L values of 1A p 
proxi- Krom (2.22) it follows that a 0) unplies Saal, —VU 
Ol ap- arars : si () 219 for all a Sut s, 0 implies a2 =0 Hence a 0 
vect unplies ad () Thus our vectors Ww, split Into two 
; a Interchanging @ and y and using (2.18). we find subsets: The subset S, consists of 1 orthonormal vee- 
tors w , ws, (ineluding w,) having only com- 
1 Tt _ . 
K pl p10 Tay pls ponents in the L-+-1 directions u, — Pha 
=. la , é- 290 subset S, consists of \4—L—1 vectors orthogonal 
E th soins to u 
2 19\! In particular andl = _— It follows that 4 mast be either 1 or L—1 In 
~ le I i i , ati ‘ HhHIsSt) ULLCSS a » ( . . 
. wi . the latter case, the 1 vectors ws in the first set will 
yor y is | ; 
ctors Letting a=y in (2.19), we find be the eigenvectors u,,, » Ui, themselves In 
vah particular, w, is uw, so that the corresponding max1- 
e pl Kak a ( eS mum of (w u uw u is zero. This occurs if and 
fviy only if « N The more interesting case is that in 
, : - which 
IV! . We substituts 2.20) and (2.21) in (2.17) to obtain — 
L=L. 2.28 
eet \ wren 
‘tors , ws Ai K p)A,Aq4, =V ax p The term in the braces mn (2.26) vanishes for 
tl 299 j ; ; j ; iy except W hen 
her ind p. This represents a set of 1 linear equations in 
the 1 unknowns (az ag? )2. It can be solved 
r explicitly to vield 
rs \ \ a a S aaa 
p A Tt va \ 
a I] I] , » 90 
28 () | N. (2.23 a: ie ia 
214 If the value i () Is compatible with the con- This result can be checked by the Lagrange inter- 
“straints (2.14) and(2.15), the minimum value of (a polation formula [17 A similar solution of a closely 
s clearly Zero We suppose for the moment that this related equation Was given by KK Loewne! Math 
snot the case. so that Z. 38. 180-181 (1934 
2.15 We now let B=8, in (2.25), multiply by a;, and 
a? <0) 24 sum $v (2.14) we have 
_ P P 
Then we can solve (2.23) with p tors In terms >» a rT) (a 0 
16 ol the a We can also eliminate s,q between (2.22 A : 
| wit! p and any othe value of to obtain 9 30 
\\ ' 
a \gain, we can solve explicitly for the 1 unknowns 
0 r 9 or . 3 : ' 
Updg3l " a, (a 2.20 a i.xp. Using (2.25) and (2.29), we find that 
7 ; ) ; / 
ints Ly x, Ett, r 
Substituting this and the value of s In (2.23) we ds I - J 2.31 
obtain | Ka, — Kp, r=t Ay—A 
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The solution (2.31) formally satisfies the conditions 


(2.14), (2.15), and (2.17). (The square roots must 
be chosen so that a B42 is (A Kg times a function 
of i, only times a function of 8, only In order to 
be admissible, the coefficients must be real. It is 


easily seen that this is the case if and only if 


K3 r K3 \ a 


] 2 l l 


The vectors ws in the subset S, are orthogonal to 


Ws, - - -» Ws, and u It follows from (2.31) that 
they are orthogonal to u,;, , Ui, This is pos- 
sible if and only if to each 7 with 717 _ t3 

there corresponds a separate xg W ith B¥B > Bs 


such that «3>dj,. 

We now consider the possible minima of (a%) 
Choosing a particular set of 8, and 7, we find from 
(2.31) that 


\ K K r 
4 8 
(a>) I] Il 2. ded) 
B,=p *p KB, i,eI d, r 
Beeause of (2.32) 
(A, Kg )( Ky r, ) (Ky A, ) Kg A; ) 
| -- <i (2.34) 
(Ky Kg (A, A ) (Kp Kg (A, r ) 


for 7, >p or B,<p. 

Thus, the right-hand side of (2.33) is increased 
by dropping any pair Ai, KB, This means that its 
minimum will be attained when the sets 7, " 
iru1, By. - - -, By are maximal with respect to the 
properties required of them. We further note that 
the right-hand side of (2.33) is an increasing function 
of d;, for 7,Ap and a decreasing function of xg, for 
5, = p. 

Keeping these facts in mind, we construct the 
.» tray, Br . . « Bz, Which minimize 
(a?)* for fixed unequal d,; as follows 


sets 2%, .. 


Let 

; l 2.0 
Let 

lo—=min! 1} ABAAY< Kg<c A 2.36) 
Let 

B,=max} Bixg< X,,}. oe | 


Then inductively, let 
tyg =min{ t|Hks 3 Ai Kar— 2 29 


B, max: 6 ke Ai, 


Because of (2.11) the set of 8, includes p If P is 


not included in the we Can easily construct 
solution of the eqs 2.14 and 2.15) with a () 
To do this we define the (2 by 2.01 The 
a. Wy Uy») 0) 

Conversely, if (2.24) is violated so that the minj. 
mum of (a2)* is zero, we can consider the problem of 
minimizing a coefficient (a with q<—p under the 
constraints (2.14), (2.15), and a 0 for some set of 
i<_p including p Sumnce kK IN _ not all the a® can 


Vanish Therefore this PELE problem W ill havea 
non-zero solution for some q and some set of i<y 
The minimizing conditions ugaln lead to the deter. 


mination of sets 7, and 8, by (2.35), (2.36), (2.37 
and (2.38) The Integer pis included in the 6, but 


not in the 2, 

Thus, a2 ~0 if and only if p is one of the 7. | 
follows from the construction of the 7 that this 
will be the case if and only if there is a «xg such that 


Np-1< ke Ay. This is assured by (2.12) for any } 
satisfying (2.10) Condition (2.12 is therefor 
necessary and sufficient for (a2)? to have a nonzero 
minimum 

The minimum value of (a2). is now given by 
(2.33) It is a continuous nondecreasing function 
of the eigenvalues \ Hence its minimum with 


respect to the XA Satisly ing 2.10) will occur for 
A.=A We may remove the assumption that the | 
A, ure unequal by a limiting process. This will alter \ 


the inequalities in (2.37) and (2.38) slightly 

As we pointed out in the introduction, we can 
replace the matrix A by a symmetric operator on a 
Hilbert space We need only assume that A has 
at least p discrete eigenvalues defined by a minimum 
maximum principle [1,5,19]. For then the first N 
eigenvalues of A and _ their corresponding  eigen- 
vectors are uniformly approximated by those of ar 


are 


l g ~ N ure lowe! bounds for the first N 
eigenvalues of A, there is an ex which goes to zero 
as N’-—-o such that A éyr are lower bounds for 
the first N eigenvalues of Ay Also, Aw—ey- IS a 
lower bound for the higher eigenvalues of Ay 
We let N’-—o@ for fixed VW and N Using (2.33 


and (2.34), we obtain the following theorem 


THEOREM 1 Let «,<Ko‘ <ky be the 
Rayle igh-Rit- upper bounds for the first M_ of the 
erge nvalues <o* of a symmetric linear 
operator A Let A, * < Ay be lowe? bounds 


for the first N edge nvalues of A. N > V/ 


Let 


Define the numbe rs 


i h, by B,<B B, by 


/ min:?7A 
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N’ matrix Ay’ for N’ sufficiently large. If \ 


the 


f pis 
uct 4 
1? —() 


Ther 


mini- 
eM of 
T th 
Set of 
> can 
laVe a 
tSp. 
leter. 


9° 


md 


3, but 


ne 

this 
| that 
ny i 
efore 
nzero 


n by 
ct ion 
with 
r for 


} 
t the 
alter } 


can 
on a 
| has 
mum 
st NV 


iven- 


of an | 


st NV 
ZeTo 
s for 
is a 
Ay. 
>) yn) 


«eo 


the 
f the 


near 


unds 


If | 


/ 


\ 
; 


B max 6 K3SA 


/ 


; min 7 ABs di, <Kg<A 


B,=mMax!B kg<ry_,, (2.40 
is the normalized Rayleigh-Ritz e1gen- 
bound Ky and Uy W the 
corre sponding to the 


Then i} Ww, 
ector corresponding to the 
ormalized ergenvector of A 


(je nvalue vA, . 


- 3 
\l—3 u u.. uy) | 
h y I h r, \ K X ) 
; l—= = T |— _—_ 
Ai, Xr, Pipe, (h h ) (A, A;_) 
(2.41) 


The right-hand side of this inequality approaches 
zero if either x, approaches A,,; or k>_; approaches X,. 
It is near one if the error x,—X, is small relative to 
the approximate spacing Ress X, and if the prod- 


ucts of errors (x,—A,)(ks,—A,,) are small relative to 


the products of approximate spacings (x, K3)(Xp r, ), 
If lower bounds X, are not given for all the eigen- 
values 4, we can always use a lower bound for any 
particular eigenvalue as a lower bound for a higher 
In particular, we can let \;=\,,; for i>p 
and d, for 4: Pp If p>|1, (2.39) requires that 
r, Then L=2, 4=1, B,=p—1, »=h=p, 
p+ This leads to the simpler bound 


u, Ww, u, | 

s 4 { (k, X,) (x ») | 
> lh > f (2.42) 
r, (ny—&y-1)(X,—da) 


Even though this bound has fewer factors than (2.41), 
it is, in general, smaller. 


eenvalue 


For p=1 we can take \,;=, fori >2. Then L=1, 
»=B,=1, 2=—2, and (2.41) reduces to (1.6). 

Example. We apply the Rayleigh-Ritz method to 
the matrix 


0 i Ze 


using trial vectors with vanishing third component 
This amounts to finding the eigenvalues and eigen- 
vectors of the 22 matrix obtained by striking out 
the third row and column of A. We obtain 


(2.44 


with the corresponding vectors 


(2.45) 


If the first two coordinate directions are replaced 








by the w, and w, directions, the matrix A becomes 
r em 
K 0) 5 Ki a 
Lowe 9 AR 
i*=| 0 K =i," (2.46) 
I l 2 ‘ 
5 K —ske'* 23 
. ¢ “ A 


We now obtain lower bounds forjthe eigenvalues 


bv means of Hadamard’s theorem on determinants 


[16] 


(2.47) 
For \, we use the fact that the largest eigenvalue 
must exceed the largest diagonal element. Then 
(2.48) 
Substituting these values in (2.41) we find that 


u,) <0.08025, 


Us) <0.09980. (2.49) 


In this problem one can, of course, determine the 
eigenvalues and the corresponding eigenvectors ex- 


plicitly. We find that 
u u, U u,)=0.000328, 
W2— Uz, W_y— Uy) = 0.002238. (2.50) 


Thus, the error bounds (2.49) are a good deal larger 
than the errors themselves. 7 

Our error bounds depend upon the lower bounds ,. 
These were chosen rather crudely and could be im- 
proved in various ways (see, for example [5]). In 
order to determine the effect of such an improvement, 
we replace the X, by the eigenvalues A, in (2.41). 
We then obtain the bounds 


u,) <0.000330, 


u») <0.002243. 


221 








These are very close to the actual values (2.50 
The simpler bound (1.6) gives 


O.1L95136 y Ov 


if the values (2.47) of A; are used, and 


lh u uw u OOP 364 , aes 


i they are replaced by the A, themselves In 
that the bound (2.41 
better than the simplified bound (1.6 


hot] 


eases we see is slienificantly 


3. Neighboring Eigenvalues 


The condition (2.39) implies that the eigenvalue A 
is simple If this is not 
eigenvector uw, is not uniquely defined 
has multiplicity m, uv, may be any 
If m 1, there will always be such aw 
orthogonal to w,, so that the minimum of wu 


the corresponding 
In fact, if X 
element of an 


the cause 


m-space 


is zero 

We must reformulate our problem. We 
PIN 
that element of the m-space which best approximates 
wv This ul 
of eigenvectors corresponding to aA 

The condition (2.39 implies not only that A Is 


value of (w,.u when wu is taken to be 


Is the projection of w, into the » “space 


and \ re 
That is, the 


separations 


simple, but that our bounds «,_;, «,, A 
vood enough to reveal its simpli ity 
error in our bounds is smaller than the 
between \,_;, A,, and XA If this is not the 
we cannot distinguish between a simple and a mul- 
tiple eigenvalue 

Suppose now that the uppel and lower bounds for 
has 2 , eae 


close together 


Cause, 


show these eigenvalues to li 


Suppose further that 
K : a \ 3 


so that A and 
cluster of eigenvalues about X 
of approximation \ like an 
multiplicity m. We ask well w, can be 
proximated by a linear combination of unit lengt! 
of the eigenvectors W,. u u 


are known to lie away Irom the 


Then, 


to oul degree 


behaves eigenvalue . 


how aup- 


This problem is equivalent to that of minimizing 


>> (a under the constraints (2.14 2.15), and 
2 10 By (2.20) we have +r, . 0 unless a=6 
oraor sp Pp \loreover.r, K aX (), () 


for a#p. This means that we would obtain the 


same minimizing conditions by imposing only the 


constraints (2.14) with a=8 or a or B p and the 
single constraint (2.15) with a=8=p. The latte 
mav even be replaced by the inequality 

=: A ad ’ > 2 


and (2.15) determine 
but thre 


The other constraints (2.14 


which local minima actually occur, local 


SeCOCK the 


minima themselves are determined by 2.14) ay 
>” 
eds 
The same situation applies in the case wm der con. 
sideration here Necessary conditions for a mani. 
mum of : > a are determined by the constraints De 
2.14), (2.10), and 2 ‘J 
Let / bye the whiit Vecrol i] thre directio} ot t} sél 
projection ot ww, mto the space spanned bw yw | 
uf 
> the 
/ S | f S ( / 
if 
let / / pe orhelr li ear COMbDNAtIONS ¢ rc 
i : i Suc Lhint al ] 
orthonormal Let { ta DP, p yy a 
and pul 
4 
Then by constructiol 
! (} PA) Db 
and 
lie 
a = ] , f A 
Cond TI 
\loreove! 
= 
| | 
ee" ad : ee / A i ee” 7 & 
Thus, if we let 
f 
\ / iY 
\ np 
\ Ol rw is > 
\ 
We have from (3.2) and ».¢) that 
i et 24 
ip 
Ssrmce the wu nre ortho ) i| 2» 14 Hpecomes Po 
\ 
he 3.10 
Qur problem is thus reduced to minimizing 
Inder the constraints (3.10 9), (2.10), and (3.8 0 
The conditions lol local 1 tiitiaik found to be as ll i] 
SectLlol However, the constraints Ri togethe 
with the facet that » \ relegate the vectors 
lh : , to thre “Tas S orthovgol | lo 
Furthermore, the conditions 5) eliminate tl 
values p 1, p 1 from the 
rr. , W 
Thus we find the follown i theorem 
THEOREM 2. Lef <ky be the Ray- 
/, igGh=-hit ippe hound hie hist \/ eidel alu 
\ aa * of ad hie miitia? operato | La 
te = Ty 
\ A. hy lowe hounds fol the first \ P [Un | | 
, t} 
een ale kon a cert pa f 9 let 
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and 


Con. 


man. 


raints 


ns of 


an elastic 


The constant € is positive and small 
+} 


h A 
k r 3.11 
Dh hie the numbers My, / —_ ‘ — = , Ba 
hy 2 $() with \ > A eliminated from the 
cet of As, and 5 h eliminated from the set 
Th 


Let w he the Rayle iGh-Rit that (GIVES 
the bound « Then the re exists a linear combination of 
unit Le nath i, OF the = _ = of A 
euch that 


eige neector 


eigen rector \ 


As in section 2 we obtain a simpler but a weaker 
nequality by putting Ni=A for + >p+m and 
=A, for 7 
This leads to 


p when p>l 


For p | we only have to put X 


obtaim 


l 
3 
By the same reasoning we can show that there is a 
combination of Up, , iu that 
with O q m We eliminate 
except for « from the x, and 
Noi-ts , A from the A, in forming the 
and Then we obtain the Inequalities 
and (3.13) with w replaced by w u by il, 
by 
EXAMPI I We consider the vibrations of a 
lorm beam which ts free atts ends and which lies on 
foundation with small linearly varying 
It satisfies the differential equation 


neal 
approximates ww 


Sets 
3.12 


unl- 


elastic constant 
y*’ eru AU O< y | 3.15) 
with the end conditions 


() 3.16 


We 


ie scalar product 


introduce 


uvdr 3 17 


and define the symmetric operator A by the bilinear 
form 

Au, v u’’v’’ +-euvlda 3.18 
If the eigenvalues A, <A, of A are defined as 
the successive minima of the Rayleigh quotient 
Av, v)/(v, v), they coincide with those of the problem 


3.15), (3.16). 

We are concerned with the two lowest eigenvalues. 
For e=0 they both vanish. Corresponding eigen- 
functions are 1 and 6°'?(1—2r). We 
trial functions in the Rayleigh-Ritz method 
obtain 


use these as 


We 


P | h € 
») 
l 
h l 6 € >. 10) 
) 
with the corresponding vectors 
u 2 | 6 - l 27 
w,=2-'/711+67-'7(1—227)|] 3.20) 
To obtain lower bounds A, we note that Ap, v) fis 
vreater for e>0 than for e=0. Thus, the A, are 
eed . oo . 
hounded below by the eigenvalues A, ot the problem 
3.15), (3.16) with e=0 These can be found ex- 
plicitly. We find 
X, =r =0 
KX. = 500.462. 3.21 


Condition (2.39) is violated so that we cannot say 
well w, approximates uw, without improving our 
However, we 2 state 


can use theorem 2 to 
are linear combinations @, and @ of u 


how 
bounds 

that there 
and w such that 


211 1—0.0005912¢ 


2) 1 1—0.0014070¢ 


Thus we have shown that w, and wm) approximate 
linear combinations of uw, and Ue in the mean square 
SCTIS¢ 

When, as in this example, A is unbounded, it 1s 
often more desirable to have a bound for the devia- 
tion CAGe iy) Ww i,) rather than(w,—a,,w,— iy») 
In order to obtain such a bound we note that the 
quadratic form (Ar, c(v,r) is positive definite for 
( r 
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Hence we can define a new scalar produc { 


la0,2 Au.v) +e(u 2 93 


A 


We how define the operator A by 


A 


Au,i u, 3.24 


a 
Then the eigenvectors of A 
A, and its eigenvalues are A Cc Applving 
2 A and expressing the results in 


theorems 1 an to 


terms of A, we find the following 
THI OREM de l nde the hypothe Sé8 Of theore Wi ? there 


exists for any constant e X, a linea combination “Woy 


u DD. * -» U; . such that 


Au.u) +eluu K ( 3.95 


When the multiplicity m=1, % is a multiple of « 


and we have the analog of theorem 1 


Applving theorem 3 to the example 3.10), (3.16 
we obtain the inequalit, 
gl ( é/ uw u dy 
D6 1. 704e—<¢ ( 1.704. | 0.00059 12. 
» «yj 
>. 6 
for any ¢€ > (), The function u depends uUpoT ¢ how- 


ever 


4. Pointwise Bounds for Eigenfunctions 


When A is a differential operator, theorems | and 
3 give bounds for the mean square deviation of the 
approximate eigenfunction w, from the exact eigen- 
function u,. It is often of interest to determine the 
value of uw, at a particular point 

In certain cases a pomntwise bound for the devia- 
“iat a point comes directly from the bound 
For example, we show that fo 


tion w, 
(3.26) of theorem 3 
anv O<¢<1 


wm, (& u(&)\** G(éé Mu Mu elu ] dy 
R 
where 
Git A Isinh?s / 
WLe.e) , “IsiInh’d Sin-o 
hers g | 


are multiples of those of 


sinh 6 cosh 6 Sin 


2 sinh & ecosl h | 2¢ 2? sil heosh | ) 
sinh h eos bh Sit bh cosh h Sil h 
2¢) sinh (1 23 sinh 6 cos b 
sit) h cosh h COS h | 2§ cost h | 2¢ 
| 
h ( ( $9 
i) 
Hee, fol example, ‘ Thus, 2d YIVeS i boun 
for |w,(¢ “(sg 

In the case of partial differential operators sucl 
bound may or may not exist. If it exists, it is diffi. 
eult to find. 

However. one can use the following adaptation of 
the method of Diaz ana Greenberg 79 cl 15.18 
For the sake of simplicity we present it only for th 
ense ot a special second orde! operatol 1 two dimen- 
SIOTIS 

Considet! the eigenvalue problem 

Aw rer.yyu AQ reyju } 
on a two-dimensional domain J with smooth bound- 
ary ¢ Here Als the usual Laplace operatol Th | 
function q Is positive iil dl is non-necative, and hot! 
are ContInUuOoUSsS In the ( losure of 1). The boundary 
consist of two parts ¢ and ¢ and we have boundar 
conditions 

i {} ol f 
Ou 
1. / / / {)} Ol ( / () $4 
On 


We alar product 


define the se 


on the linear vector of funetions which ar 


plecewise continuously differentiable In Dand Vans! 


space 


on ¢ 
Let u 
ine to the eigenvalue X 


be the normalized eigenfunction correspond- 
and let w, be the funetior 
corresponding Ravleigh-Ritz upper boun 
h for A; Theorem | YIiVeS il bound lor the deviation 
In Norm (wv “Uy, W Y 

We wish to upproximat the value of w, at a 
interior point of D. which we choose as the origin of 
oul system We the fact that 4, 
satisfies the differential equation (4.3) with A=A, 
the boundary conditions +.4 Let T(r. y) be a 
parametrix for the differential equation a satisfy- 
1.4 That Is, 


to a 


coordinate use 


ana 


Ing 
4.6 


where 9? is any Lwice continuously differentiable 


function such that 
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ia (} on ( 
2¢ or 
’ LT —0 on f 1.7 
On ; 
Multiplving (4.3) by T and integrating by parts, 
"( find that 
19 ; "u o. 7] Ad } A q P \dady 1S 
OUT . ; 
ind Replacing u by ul and X by Kk, on the right, we 
obtain the value 
uch g 
| diffi. di 
wm (O00 1h Ad / K pf] Cidardy 19 
iOn of “D 
18 
mr thi which cun be computed by quadratures Using 
imen- | Sehwarz’s Inequality, the normalization of a and 
he triangle inequality, we find 
) 0.0 u.(O.0)) <x 31 y 
} 
il i uw u 
yund- 
The ( } 
bot} z Ag / Kp) dedy ¢ £.10 
iry ( “D 
dar 
Thus. the bound (2.41) for (wz Un», W u, together 
with the bound «,—A, for x,—A, provides explicit 


ipper and lower bounds for u,(0. 0 These bounds 


‘ 


lie close togethet if the error bounds ? 4] and Ky nN 
14 re small 
The same method applies to the function @, of 


theorem 2 If 


: y S* alu 1.11 
15 mand 
we find from (4.8) that 
ure 
LHS — 
/ o. 0 Sa uf Ad j Nd C\drdy 
ond- ~ PD 
pone 1.12 
und 
tion | Hence if we again define the approximate value w by 
19). we obtain 
aul 
n of ; mg 
/ ()} {) M {) (} WwNX A X i K 
Us 
and 
ye | | Lh / Li / 
sfv- 
: °— ) 
q Ad / K yf ty dxdy ¢ . . 13) 
1 6 “_ 
This Inequality together with theorem 2 vIVeSs 
able 


ipper and lower bounds for a »(Q0, O 


J. B. Diaz 
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Sequence Transformations Based on 
Tchebychetf Approximations 


John R. Rice ' 
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LOo0 


ireal sequence and define the as follows 


ransformed sequence C P 


Consider the segment P, Pr Py Po of P and determine a,,, b,, and c, so that 
IAN m_? c, P, n=m, l, 2,7 3 is minimized The mth term of ¢ (P 
ike be The {Tec of the transformatior C sequences of the type 
p > Ja \ p ] > aon S hon’. and Por s considered In each case C is 
( o be very effeetiy i accelerating convergence or decelerating divergence For 
imple if the second sequence behaves as an is > then the transformed sequence 
behaves as ar \ similar transformation ‘D is defined by approximating in the Tcheby- 
ff sense a segment of P by ab, cos (6 no ( The effect of “D is studied for sequences 
f tl ibove type and also for p >“ ,a,NI cos (0 These sequence transformations 
ire similar in nature to Aitken’s 6-process and its gene ralization \ comparison of the 
two types of transformations Is rasa Several « xample s are given to illustrate the effect 
of € and ‘D on various sequences 
l. Introduction lf b 1, then Lim /(A,,,”) is just ¢,,; if 5, | then 


this limit does not exist, but the value ¢,, is still taken, 
c,, being the antilimit. 

Explicit formulas fore inh terms of the p have been 
) obtained in 5] These formulas are not linear in the 
p, and the transformations exhibit some marked dif- 
ferences from the classical linear transformations. In 


Let P Pri Pe ..> be a real sequence and 
let T bea sequence transformation with 


O=Fi(r t,( Py,p » Pa)|? = 


\ large Variety of such transformations have been 


studied, many in the modified form of summability 
ethods for infinite SeCTLes, which are designed to ac- 
elerate Ol induce the convergence of the sequence r 
Such transformations Trea be based on upproxitat- 
ya segment of 2? by au function of a given form \s- 
that p behaves approximately as f(A.) where 


| represents 


i! " 
SUTTIE 


determined 
2, The 
Piety be determined so as to minimize 


Sore parameters to be 


Consider a seement p i.2 7 of P? 
mrameters | 
e errol of the approximation Ink SOMME Sehnse 


A,,,”) 1s 


swice 


a known function of m we may set 


Pr, ip Lim /(A 

The transformations of this paper result from as 
SUMING p to behave 
The transformations are denoted by C and ‘D. re- 
spectively The upproxination to vive the mtl 
erin of C P) and ‘Di (2?) is based on four or six, respec- 


aus ab Col ab 


COs (0 iQ ( 


tively, terms from 2? The parameters are determined 
sous LO Tinie 

max a h ( ) } 7 Ti 7 2 9 ) 
OT 
Max d h COs( FA lio ( DD | ¢ } Wen eo 0g o 

| 1 ( W \1 

{ 
\ 


many cases C and ‘D are much more powerful and in 
other C2uses r Thiet converge quite rapidly and the 
transformed sequence Tay be identically infinite \l- 
though C and ‘D are not regular they do have a prop- 
namely that if both ? and 
the transformed sequences converge then they have 
the same limit 

It is not surprising that C and ‘D are 
ces of exponential type 


erty of }0 nt CONVETGENCE, 


effective for 


sequel 


‘ 
_ a r fi . . a ° ° | 


lf A \ then (CU?) converges or diverges as 
\ instead of rN ‘D(P) behaves as (A ol 
As/A ‘D is effective for sequences of the more 
veneral type 
) ) 


a Pot Ja \ 


COS “7 /L@ fi 


and Cis useful for such sequences In some Cases 
These transformations are also very powerful for 


alter hating rational sequences 


“¢ 
4 Pp | "> 2 n > 5b n'n > 2 . i) 
\ 
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It is shown that € accelerates the convergence by a 
factor of 1/n That D accelerates the convergence 
even more can be seen from the examples although 
the actual factor of acceleration has not been deter- 
mined. An example shows that ( and “Dare very 
effective for a sequence of the form 


yg Z, } 


The basic principle of these transformations 1s the 
same as that of Aitken’s 6°-process and its generali- 
zations. € the direct of the 
Shanks [S| has given the most complete account of 
this theory. € and “D have many properties in com- 
mon with these transformations. Tl 
and “Din Shanks’ work have the 
and “D for sequences of type (1 
Cand ‘Da 
6 -process accelerates the Conve! 
gence ot (3 by a factor of only 17 C and ‘D are 
much more effective for the sequence (4 There are 
sequences for which the 6 -process ana its gvenhe rali 
a general 


analog o--process 


Is 


sire 


powe! 


‘) 


ana és lor Ssc- 


quences of type » re definitely more pow- 


erful Aitken’s 


zations are more effective, but as 
and “‘D are as effective or more effective than 1 
analogs in Shanks’ paper. 

A comparison of this paper with Shanks’ will indi 
cate that C and “D have many properties which have 
The purpose here hh; 


neu 


not been developed Prheain is 
been to establish these transformations as tools of 
numerical analysis and not ol analysis When read 
ing Shanks’ paper one has the feeling that much 


remains to be done in the study of these transforma- 
tions and that a complete understanding of their 
behav lor for sequences of real and comple x functions 
will require a penetrating analysis 

It is typical but somewhat disappointing 
transformations may effective 
problems and least effective for difficult problems 
If one is solving three simultaneous linear equations 


that thes 


he most ior ensy 


by Gauss-Seidel iteration and the characteristic 
values determining the rate of convergence ure Qh 
5, and .1, then these transformations increase the 


rate of convergence dramatically Sut if one has 
LOO equations with 10 characteristic values between 
99 and .995, then the elimination of even the five 
largest ones is not very On the 
hand if the largest characteristic value is .995 and the 


others are, Say, less than .9 the transformations are 


significant ovhel 


very effective in accelerating the convergence 


2. Derivation of the Transformations 
2.1. Preliminaries 


sequences are denoted by P. (J, and their ele- 


ments by Pns ¥ A sequence transformatior 


is denoted by a script letter as C and “D Braces 
denote a sequence or set and (4g Is rea 
“the set of « such that We denote by Ola 
re VN _ 

a Is S Ix << Is S 


The analogs of C 
as C 


cuide C 


and o ] two functions suc! that 


and, for x sufficiently lare 


() 
constant 
Let P be aw sequence ol real numbers \ trans 
formed sequence (J Ini? > . ot P. is oh 
tained as follows: Take segment, p m, 
nN i of P?, and upproximat the values p is 
function ol hy ab C Ol ah Cost Ho c TI 
corresponding ¢ value Ss then assiened lo g 
Kor these particul siormations to be effeetive 
Lie sequence P | lst } SOT SeTs¢ behave (X- 
ponentially \s i! CNXal rl ol Sue sequences 
consider the rt il J tors «ce ned hy 
| 
whet Lisareal / Let be the solutio 
ol 
| / () 
Phen | The sequence OF vectors 


ned will converge to * i A tends to zero 


‘nds to infinity lf Lis a normal niatrix wit 


us wv te 


Characteristic vectors ina alistinet real Character- 
stic values A,, ther 
| Ja \ / 
If the Ay are compley ther 
| S y \ Cos 4 ren) / 


2.2. Tchebycheff Approximations 


We wish to approximat segment of P? by ab"+-« 
The ise of least squares approximation appears to 
be unpossible due to the difficulty of the nonlinear 
equations involved ly y} thre theory ol upproxima- 
tion by ab cin the Tehebyeheff sense is cl veloped 
in detail. The follow result is give! Let p 
p ,~p and p be four consecutive values of P 
to be upproximated Then ¢,, of the best aupproxima- 
tion. da h C,. IS Give Dy 

/ p }} p p “ 
f ) 
2(p p / p 

Let p... p be consecutive values from P 

to be approximated by ab"cos(@+ n@ d. The theory 


dis 
of th 


r@ 


d 


nation hy ab Cost 
fol 


ol Tchebve hefl aupproxin 


not complete and uh ¢ formula 


x pli il 


best upproxination has not been rigorously estab- 

lished. However it is conjectured that the following 

formula is valid. Set p p 1=0, 1, 
t+. Then ) 

_ ‘ R : { s Q 0 

- | 


10] 


‘ 
A 


oe) 


V1} 


This formula is based on the assumption that the 
hest Tchebycheff approximation Is characterized by 
five times of the error function. 
see 6] lor a discussion of the characterization of 
hest Tcheby cheff approximations 
Note that (6) may be written in a much simpler 
form Let s ‘ %), @@=0,1,3,4; then 


ie alternation 


\ \ ) 
lrans- 
Ss ob- 4 ‘ > 
2.3. The Transformations ( and ‘D 
iS | 
Thi Let P? be given, then a new sequence @ 1s found as 
follows: Kor each m such that the denominator of (5 
etive s nonzero g c,,. if the denominator of (5) is zero 
eC eX und the numerator does not vanish then we ussich 
ences if both the numerator and denominator are 
ero ~P This transformation of P? into J may 
e written in operator form as 
() F Pp 
utiory : ( " 
A transformation ‘D is defined in a similar manne 
vy (6 
‘CLOrs The repeated application / times of i and ‘D is 
» zero! denoted by C*(2’) and ‘D .? 
Wit! 
ic tler- ; , 
3. Analysis of the Transformations 
3.1. Algebraic Properties of ( 
The study of the properties of these transforma- 
tions begins with a simple algebraic property of C 
If P p thenaP bis defined as ap h TL 1? 
The following result mav be established by 
rect computation from (09 
)" c : . = 
rs to lH OREM |] ( ap h af pP /, 
neal 
ima-} Itis not true that C(?;+P co C(P 
oped 
> D A . 
. p 3.2. Transformation of Exponential Sequences by ( 
na- 
Let P be a real sequence of the form 
7 {> p) T >) aX; i 1.2. } 
a P 
eory vith \ A If A e P does not converge but 
ad liverqes froin p In such n case p Is the antilimit 
of P. This sequence Is of an exponential type and 
ordi ( should be effeetive in increasing the rate of con 
wing 


vergence Ol decreasing the rate of divergence, 


0,1, Let P" Po; then 


C(P)=C(P)) +p 


The general term of C(/” 


puted to be 


may be explicitly com- 


>) ajar N,/A1)™ (1 +A,) (1+A,) (A A,)° 


(A;/A,)" (1 —A,—AI+A 


The denominator of this expression may have only a 


finite number of zeros as a function of n. We may 
write the denominator as 
a, (1—d, —AF+A3) +0(1 
since A,/A | for 7>1. 
The largest term in the numerator is X° Others 


among the larger terms are 


Thus C is seen to eliminate the largest exponential 
term from P’ 

The repeated application of C will eliminate the 
largest remaining term. The largest term in (7(2” 
= 


Phe above analysis has established: 


Tueorem 2: Let P={pot+DdJa,r’ \, \ 
the) 
( ‘ 
Nu > ld u,)’ M | 
CC Por 1+0(] ui<1, i>. 
\ 


Ss See) 


Thus € to have a desirable effect on se- 
quences of this Lype, which was to be expected, In 
the next that C is effective for 
some sequences of a completely different nature, 


section it 


Is seen 


3.3. Transformations of Alternating Rational 
Sequences by ( 


Let pP be a sequence of the form 


{» Por 1)" Slam Stan}. 


We shall consider € (?). 
First take the special cuse 
"an 


r p | 


The ath term of C(P) is 
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After some manipulation this term is seen to be of 
form 
ah 
Port 


‘s 


nh 


Porta 1)", k : Of )| 


Now consider the ceneral Situation p may 
written as 
p p a,/ T Ja fi j j j 


for 7 sufficiently laree It is clear that the den 
nator of ¢, is 8 L)"kn O(n since it is In 
in the terms of ? The denominator of « has ¢ 
a finite number of zeros as a function of m 


terms in the numerator are of the form > a.aa(X-r-,)) 
ajaln fi l fi 2 / ) COS a Q ps da rl H 
; and 
aja;|n'" n | } a n } 
+ Ber MAAR "G 
2aia,| i 1)’ i ar i l / 2 
By direct calculation it is seen that this term is where H,;, and G,, are functions of the 2’s, 6's, 
On? ‘) @s that are independent of ? /] and G 
Thus we have established trigonometric functions with coefficients depending 
oh powers ol the \’s 
THueorem 3: Let P p 1)*> ** an't/> <*. bmn: An examination of the above expression and 
hen fom om sufficiently large ae: “aa fam ol C simplification of notation leads to the following 
P) is theorem 
THeorem 6. Let P >i 1a ,X" cos (0 d 
Pot >a " \ \ Then t) h term of DIP 
where k=k,—k2, a=—a,.k(k—1)/86 MS >) bu? cos (a 
Examples of sequences which are approximately ol p 
this t\ pe are c1Ven bv the partial sums ol alternating as 
infinite series of rational furctions of 7 Qne car 
consider such sequences directly and establisl tin) 
analog of theorem 3. eCLce pt fOr poss hlyaft her of terms 
Qne ma nlso establ t! following result 
3.4. Regularity and Joint Convergence 
A transformation T is said to be regular if. the THEOREM 7 Li [? >a Ni Sy \ \ 
existence of Lim p, implies Lim ¢ Lim p, wher FE om the oth, & Di} 
G=F (ir Since € is a transformation designed to 
accelerate or induce convergence one would like lor 
C to be regular That - is not regular is seen by A >) 
the following example 
\ max (Ag, A5/A 
P=i{ 9, |p l+1/k, p l+-1h 
ercepl possibly af A ff, 
Every other term of C Is Infinite vet P cor No attempt has bee made to analyze the eff 
verges. Although Cis not a regular transformation of “D on alternating rational sequences: however th 
it does have the property of joint convergence \ examples indicate tl D is quite effective for 
transformation ‘ is said to have the property ol creasing the rate of convercenes 
jornt convergence if the convergence ol r and T r It Is simpl ieaatt to construet sequences 10 
imply that ? and T(/?) have the same limit Thi which ‘D is not recula tis probable that “D hi 
proof of the following theorem will be vIven in sectiol the joint convergence property. but this has not be 


5.1 theorem 10 


ScCe 


THEOREM 4 


the 


3.5. 


Properties of the Transformation ‘D 


The following theorem Gives aleebrate 
of ‘D similar those of ( 
evaluation 

THEeoreEM 5: ‘D(aP+h 


( ‘onsidet the real Sequence 


[P? { 0 \ COS ‘iu H@ 


to The prool is by dire 


be 
aD(P h 


yTNil- 
eal with \ \ Expression t) may bye direct 
aly evaluated as the ratio of the following two terms: 


The 


established 
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the 


properties | 


ot 


} 


then 4, Comparison With Aitken's 5°-process and | Likewise, a straightforward calculation shows 
its Generalization . 

‘D(P)=&(S 

The transformations C and ‘D are similar to a T 


mily of sequence transformations first studied by reason for this relationship is simple. &, has 


been found by solving eq (S for ¢: 


\itken 1,2 These transformations have been stud 
eT Lies dl sirce by Lubkin | Samuelson 7 . and Shanks 
dire, X The most reneral and complete treatment Is Sab c i m n k. - mn i. (R 


vy Shanks A brief outline of these transformations 


I] be rivell 


Let P be given and set The fact that best Tchebveheff approximations are 


characterized by the alternation &+-1 times of the 
error function leads to the following equations: 


Ap p p 
Krom , new sequences (J Tks / k. j a >" h ( p | ( Mi HN, . i Ik t - 
rectly re formed hy 
y yp p where e is the error of the best approximation, If 
: one averages the mth and (m | )st equations, then 
as MPs oP a(i+b 
H 24 5 O7 +e=5 (Pat P 
Ap Ap Ap “ “ 
Mark, . , n+2k+1. 
These equations are of the same form as (8) except P 
s Al has been replaced by S 
are | Ay Ap : Aitken considers the application of &, to sequences 
ndir ae | | | ol the form r p > »- a;X' and establishes: 
= Ap Ap,, THrorem 8 (Aitken): Let P=) py +> Sam 
wing \ \ : the) 
Ap Ap Ap 
no aT), 7 
aad i rv l 
| Io! ui<1 i>. 
Ap Ap \ . 
This defines a set of sequence transformations &;, Iyy Thus for sequences of this type & and € are equally 
effective in order of magnitude. The coefficients of 
0) £,(P). A for &, and € are, respectively 
' Q(X \ ajao(1+-A.) (A \ 
Shanks studies in some detail the transformations ’ 
c ck © 1—A,)° 1—A 
C..C; PisGi2.4 : ] ‘ and 
é 4 k=0,1,2, qin denotes the nthtermof | ¢ appears to be slightly superior if A. >0. Examples 
al qin is given explicitly by indicate that either &, or C may be better but in any 
case the difference in the rates of convergence is 
s ld id Pn. small, 
‘ p p 2p For alternating rational sequences the two trans- 
formations are no longer similar in effect. Shanks 
hich defines Aitken’s 6 process has established 


These transformations are derived heuristt ally us PHroreM 9 (Shanks Let 


lows (riven a segment of . p , , Pp . 
effect | Ssume that this segment is exactly of the form pP {> 1)"SJam' | > Jb } 
r the da h ( Then define g,,, to be this « The basic 
pr in 7 a,. #0 h. <Oand i] } i Then the nth term o 
rinciple is the same as for ( and ‘D, but the 
, etermination of the approximation is different. 
) has | If S PL. Lin p hen it is seen that | 


O(n’ 
C(P) ES 
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with possibly a finite number of ¢ rce ptions are made with &, and é here The next set 9 


Thus €; aecelerates the convergence by un taector examples shows the effect of increasing thre larges 


of L/n whereas C accelerates t bv a factor of | A in exponential sequences The following [WO set 


The comparative effectiveness of & and ‘D may show the effects of lower order terms and the loss n 
be judged from the examples given lates significant figures Detailed comments are. giy, 
Theorem 10 shows that &) has the property ol with each set of eXamples 
joint convergence } All ot these example s were computed USING doub| v 
TuroreM 10 (Lubkin): // both P and &€(2P?)=—@ | precision arithmetic with about 20 significant decim 
converde, then digits , ot 
\ 
im p Lim ¢ . ; ~ 
: f ; 5.1. Effect of ( and ‘D on Exponential Sequences 
, The following five examples are considered 
If P converges then so does Ss p p and 
r 2 1.2 4 \ 
Lim p,»=Lim 
n n r QS Q5 5 2 
P QS Q5 
Thus theorem 4 follows directly from theorem 10 
P, y s COS / 
5. Examples 4) " COS | ay, 5 COS PAST 7 ( 


‘a and ‘D fo 


Some examples are given which illustrate various | Table 1 gives results of ¢, €?*, C’, 
characteristics of C and “D In some cases com- 33 terms of 2? P, diverges quite rapidly and ¢ Pr 
parisons are made with &; and & for more diverges at a much slower rate The factors 
general comparisons see '4] where &, is compared CP larger than .5 are 1.2, .9, .8, .72, and .54 
with several classical linear sequence transforma- P? converges as can be seen from the list Ol | 
tions for a large number of examples factors in C # larger than .5 They are .v, 6 

The first two sets of examples illustrate the general ic. So. &. and 54 The four largest lnectors o 
behavior of C( and ‘D for various exponential and | (, (72) are .S8, .72, .71, and .68 are the largest facto 
alternating rational sequences Some comparisons | of (*(7?,) is .72 
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Set of 
urges 
O Sé@lIg 
loss 0 


“iver 


loubl 
PCImM 


nces 


Drs] 
| 4 
list 0 
y 


} 


the first 


creases 


It seen that 
rratically as k 
metion of 4 Tor small» 


Is values of C r behave 
“as very erratic 


This is the ty pical behay ior 


Is a 


ind indicates that care must be exer ised in Using ng series 
repeated applications ol "i 
The largest factor in ‘D(?’,) is 9 and henee ‘DP 


onverges cs a. i 
Onlv selected values from the othe examples ure 


nven In table 2 it is seen that Cis not verv effec- ; — | 
g : * ie : ieee (j 
ive mm mecreasing the rate of convergence of / hunt 
Prarie 2 = on })'9 
/ Cu Cu Dil These sequences are denot 
respectively 
in ss te a 


sequences and tables 6 and 
and &, to be as predicted 
DA ‘D(A,), &(A 


"| is 1s due to the closeness of the two largest factors, 
more 


the next tuble it is seen that even with the lower 


and J 
rapidly and apparently ‘D accelerates the 
convergence better than re by a factor in 1/n 


5.2. Effect of (and ‘D on Alternating Sequences 


( ‘onside the sequen es of partial sums of the follow - 


951472, 


LOUNGE. 


by A,, As, and A 


ed 


The first two eXamples are alternating rational 


7 show the results for C 
by theorems 3 and 9 


€,(Az) much 


converge 


rder term deleted the convergence is quite slow 
¢ *| , 
\ Ol / _ 
en i ( PARLE 6 
PARLE 3 
/ ‘ ; ( i é { DA é i 
Cl CP ( Cl D PEE E | twit 141738 3. 1423423 
itt 1013 133 41556 $. 1413919 
{eit 121 452 11604 $. 1416627 
Sad Std RUFQ3N 13 4. 1BUON 1415882 , 14156%4 
SSS s Wi24 Js Hn4 s 41 7260 142 ) 1415946 116005 
“4 04 \" {15163 1408813 H15854 
* _ j 142 5 
~ * VAs 
The next table shows that “Dts fairly effective on 118 
p lt is to be noted however that convergence is ho 
onger monotonic even in the large range considered ’ 4 ” ° alalalalt m_ 
The denominator has a zero between the 67th and 
Sth terms of ‘D ra This same phenomenon is 
resent in ‘D P, PARLE 7 
PARLE 4 
{ COA é i Di A é i 
HO444444 BIRT 69333333 
/ cil Dil De] HO2ZTOR3 HO047619 60313724 60308043 
x 64333333 04444444 643 1 5068 64316040 
— 6930 60242424 HO314 HU313724 
{ xu ay | Rg 3 HO31972 HUS5SS074 H4U3147S U315212 
{ 608 $072 0068 61686064 60311 69285712 69314682 69314448 
‘ ; MIS s 12 Hustle HUSS4 732 69314734 64314872 
j LI “ WK)4S “ 142370 ) 4K 69300331 Ou Ti 04314621 
1563-0 693 6932539 69314 S14778 
} wud O00 j 14 314 Hus HY. { HU31 4 
“ 1 lf H5440 64315149 Pall 41471 6031474 
n " 6931439 HU3096870 60314 HY3 1460S 
133 H93 1406 IQ1STUR j | 14 
R70 1314520 60311377 { 
a - : 4 31486 69317488 
Phe final table shows that “D is very effective for 
) . . (H2RTIN4 2394 
, while ¢ Is not 169538 
Tarte 5 
I Y wu LLL 
/ Cu Diu De] 
It is remarkable that 7 can be found so accurately 
, : from a few terms of the simple sequence A;. Shanks 
ni considers €$(A,) and 7 is given to 8 places by 
~ shy al x I the first term of &4 ‘A It would require aup- 
' proximately 10.000 000) terms of A to give the 
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aceuracy 








\ll 


mre 


The third example is a ‘mixed 
of the transformed sequences diverge, but they 
pp p20, 
quences for /n 3 In this exampl C and ‘D show 

very marked lmprovement & and & It 
uppears that round off has begun to affect “D 
Although “D(A 


Seq UCHICE 


SEM ICONTET CE nt or asym ptoti Z. 586] s('- 


il oOvel 


A 


Wiis computed using double pre- 
cision, Al, itself was computed to only eight signif 
icant figures 

| E 8 
! ( é 1 D é 
The sequence | mas be modified by repeating 


each term in the sequence twice. Call the resulting 


sequence A The application ol iA ‘D, € ana 
€ to A; give some unusual results It is seen 
that 

E(A A 

C(A a, a / 2 
€.(A is of the same {ype as A witl eCucl term 
repeated twice The odd terms of ‘D(A nugree with 
those of &. ye but the even terms do hot The 
even terms of €(.4;) converge much more rapidly 
than the edd ones, by a factor of approximately 
lin. Both &C4;) and “DCAS) converge more slowly 


A,) and ‘D(A 


than & 


5.3. Three Related Sequences 


The 


quest Ons 


purpose is) to (1 following 


rate tiie 


LCs 


LW 


Isitadvar reous to rcrease the larges 


term in an exponential sequence befor: applying 
eT Is if ndvantagweous to keep the larger lerme 
away irom l at the expens ol inereasing they 
The answer to the first question appears to be some. 
times and the answer to the second uppears to fh 
ho 
The limit or antilimit of each of thes SECC UeHCES 
Is zero The first conve reves) quilt slowly and ti 
other two diverge (see table 9 ‘ 
( 
' P;, P P? C 
} 
Ps, C Ps), r ana C all converg ni 
slowly There Is an lmiprovement by un o6inctor of ( 
approximately 2 in the rane considered for ea st 
further application of C ' 
Theorem 2 states that ((7’?,) and ((2,) converg 
at the same rate while ((7’;) diverges. This 
borne out by tuble 10 ( and C produce entirely } 
different effeets for enc! SEC UCC! (C a CO 
verges Very rapidly lo ero) and te e COnVCTGES ) 1] 
evel more rapidly Errors due to the loss of . 


nificant figures have affected the 24th term of} % 
C r and the 14th term ol C P? All cal \? 
lations here were made w h 20 decimal digits. ) 
In P;, COP), CROP nd (°(?:) the first term 
the closest to zero However it is apparent that tl 
sequence formed Dy thie tl terms ol C P. 
j (), |, converges quil rapidly to thr antilin 1) 
of P? 
' 


5.4. Effect of Lower-Order Terms 0 


_ :' To illustrate the eff of lower order terms 
Che following three sequences ure considered Po l ” 7 a de ; 
considel 

P QS Q5 

pP 1.96 Q5 r = 1.2 

P 1.96 19 : 

Wiel rea \ bye compared . pP nsection ob 2 
Tat l 10) ff rf } / sS oa (yf ‘) (yf ; 
( ( ( ( ( ( ( ( ( 
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r ty Parr Eff ( é 

larges 
plyir ( / (i/ ( Cel ee C4f (*u 
terme 

then he oe 2 7 4 S40) 1 ~ “ ~ “ 

“uD D 4 % NN H, NOE . O73 j ool wt 

SOM. unned 1 On 4 

to bh 

LCTICes 


In tuble Tl it is seen that the lower order terms thre hntlmerator, the actual value of the numerator 
Meu lisappeal rather rapidly in 2?,. but the 27th terms ef | and the number of significant figures lost from P 


(P and C f’ only nvree to two digits Further 


( P and ( r are completely unrelated to € rs Tarte 12 Lo ol gnificant fig 
ind C oP Thus it Is seen that these lower ordet 
terms are significant even for € ??) and they affect Prat Size of Signif- 
he rate of convergence Of C r : I - drastically , ae ' 3 ; . ~ ; was P ~~ 
nd adversely , 
5.5. Loss of Significant Figures +08) 10 7 3 vou ie - 
oy ‘pee Bm i es ti he 9 
It has been seen that terms with a small effect on oi Hera D | ! ! 
} P can have i laree effect on C P A related 
iverg phenomenon is the loss of significant figures in the 
tor of omputation of CP) and ‘D(P’ If P is the con- 6. References 
© OU stant sequence Pn Pp ls, then C P) and “D(P 
mav be written S\ mbolically us I] A. 4 \itken, On Bernoulli's numerical solution of alge- 
iverg : braic equations, Proe. Roy. Soc. Edinburgh 46, 289 
. 106 
his , 1-1 21 A. ( \itken, Studies in practical mathematies, I], Proc 
itirely } | | Roy. Soc. Edinburgh 57, 269 (1937 
con- s} WK. Knopp, Theory and application of infinite series 
ee Hence slight inaccuracies in_ the computation ol i} Ss ar ga gy Hs ny Bo series, J. Re 
yt SIR. ( P) and ‘D aa may distort the results badly if 7? search BS 48, 228 (1952 
rm Ol S slowly Varving Qiten it is necessary to use double 5) J. R. Rice, Tehebycheff approximations by ab c. J. Soc 


Ind \ppl MIath to be published 


eal i- , . ° ra , 
recision for the computation of “D(??) even though a ia 
’ | 6) J. R. Riee, The characterization of best nonlinear Tcheby- 


‘is only given to six or eight digits 


, woah . E cheff approximations, Trans. Am. Math. Sov to be 

rm Is i fhe follow hy table shows where the loss of signifl- published 

alu int figures can occur. The first column of table 12 | [7] P. A. Samuelson, A convergent iterative process, J. Math 
r Phys. 24, 131 (1945 


efines the sequence a the second column VIVES the 


tilin “i a P Pp —- a dian s 1). Shanks, Nonlinear transformations ol divergent and 
pproXimate size of the terms ol In the range that slowly convergent sequences, J. Math. Phys. 34, 1 
he transformation is applied and the third column 1955). 
vfIVeS the transformation used The next mul 
olumns give n, the size of the terms appearing in (Paper 64B4—3s8) 
ns V 
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‘Numerical Solution of the Frequency Equations 


for the 


| Flexural Vibration of Cylindrical Rods 
W. E. Tefft 
Ju 28 LOGO 
} \ numerical solution of Pickett’s frequency equations for the flexural vibration of 
evlindrical rods, based on the three-dimensional differential equations of elasticitv, has 
| been obtained on the IBM 704 computer The ilts are presented in the form of tables 
of correction tactors to. be ipplied to the thin rod ipproximation for the fundamental 
fle ral resonance Treque ev and the first two overtones These results provide an accurate 
1" ins of ealeulati v Your i, modul is from the density dime nsions, and resonance tre- 
quencies of evlindrical rods having diameter-to-lengt ratios as high as 0.6 


1. Introduction 


Improved experimental techniques recently have 
made it possible to determine the elastic properties 
VI materials to a high degree of precision by ly hamic 
methods which have limited 
Convenient experimental 


resonance Factors 


wcurate Measurements are 


echniques are not well developed theoretically, and 
onversely, techniques based on well analy Zed 
neories are experimentally difficult A recent ex 
erimental investigation [1] ' showed that Pickett’'s 


) 2] approximate solution of the flexure problem, 


mused on Croen's 3] solution ol ‘Timoshenko’s } | 
quations may be in error by a few percent. In a 
iter paper, Pickett 5} vave ah auecurate set of 


lrequeney equations based on the three-dimensional 
ferential equations ol elasticity The purpose ol 
e present papel present the results 
merical solution of these more accurate equations 

These numerical have put the 
eOry of flexural vibration of evlinders into a more 

isable form elimination of theoretical errors, for 
s particular problem has thus been achieved 


is to 


ol a 


solutions accurate 


2. Background 


According to the elementary, one-dimensional! 
eoryv, Young’s modulus, #, may be computed from 
( flexural resonance frequency, iB of il evlinder of 


ensity p, length /, and diameter d, from the equa- 


tion [3] 
642-p/'} , 
mid 
here 7, is the resonance Irequency of the n™ mode 


| vibration and the m, are solutions of the equation 


COS Ii cosh Mi | 


The values of m for the first three modes of vibration 


re [6] 


n 1730041. 3 
iT 7.853205. | 
dl } It) QO500S8 2] 
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Timoshenko |4] added shear and rotatory inertia 
corrections to the differential equation, and Goens [3] 


reduced the resulting frequency equation to the form 


6 T pl*f: T 
mid 
/ sunetion OL dil and Poisson’s ratio, 


approaches one 


where 
u.T 
all values of Hi and lu 

Pickett [2] made numerical computations of the 
/’,, and gave approximate equations for this com- 
putation These values are not accurate 
desirable, in view of the experimental precision now 
attainable in) dynamic measurements. 
The basis for the more accurate values presented 
in this paper was given by Pickett [5] in his solution 
based upon the three-dimensional differential equa- 
tions ol elasticity 


Is i 


iis dil approaches Zero, for 


as as Is 


resonance 


\ brief verbal sketch of Pickett’s derivation ay 
be helpful in understanding the meaning of some of 
the following equations. He begins by noting that 
there three fundamental 
waves In a evlinder. Each is characterized by the 
angular lrequency p and the wavelength 27ra/y, 
where ais the radius of the evlinder These three 
combined so that their sum satisfies 
the boundary conditions on the curved surface of the 
evilindes This vives a relationship between p and 
+ expressed by Pickett’s [5] equation 12 
Value ot p there are two values of ¥, denoted by y 
and which satisfy this equation, One of these, 
‘, corresponds to a sum which represents a traveling 


ure solutions for elastic 


solutions can be 


Kor each 


wave The other, +. corresponds to a sum in which 
the “wavelength” is an imaginary number, i.e., the 
displacement is an exponential, rather than a sinu- 
soidal function of distance along the rod Pickett 
combines both of these sums to get a solution for 
displacements which approximately satisfies the 


boundary conditions on the flat ends of the cylinder. 
These boundary conditions lead to his equation aks 


the Irequency condition The present paper is 
concerned with the simultaneous numerical solution 
of thes equations 








3. Theory and Calculations 


form ol 


Pickett’s 1D results were reduced to the 
equations (7) through (12) below, by a method similar 
to that used yy Hudson {7 
an expression of the requirement that the three radial 
components of stress be zero at the the 
This equation was reduced to the form of 


evlinder. 
eq (7 by the aid of definition (10 a suitable 


manipulation of the resulting determinant to reduce 
it from third order to second ordet 


Pickett’s equation 12 


Is 
surtuce olf 


and 


2heWih) + (Ke —y*) |2W( A+ AK 2V7(3Y( A) +A 
2h kK? +4 Ke —y) WK Kk? +4 


where h and K are defined by 


where 
rv Poisson’s ratio, 
i’ Young’s modulus, 
p—density, 
a d a. d bene diameter of evilinder, 
p=2r/=angular frequency of flexural vibration, 
y parameter determined from the solution. of 
these equations, 
Yy—a combination of Bessel 
kind, 


funetions of the first 


orders ZeCTO snd One, ai cl Is cle fined by 


<~ 
— 
ty 


If we define the function 
Kk *4 : ms 4 
Sie - \ h 
hO=, A |v A ( 2) 


LWA 


la a 
K Corl ~ | 
\ i <1 
+— (4 K ] ir 1] 
~ Lan ‘ 
\ <1 J 
where the uppel terms ih braces relel lo the odd 


numbered and the lowe! terms to the even numbered 


modes of vibration. then equation 21 of 15! becomes 


ermining the flexural resonhane 
For il 


The problem ot de 
frequencies is now solved in principle 


“lve 
leneth, diameter, You 1's modulus, and Poisso) g 
ratto, one must solve 7) and |? simultaneously {0 
p and incidentally lol and ; However, thes 
are transcendental equations and involve the fre 
queney in such i compl ¥ wav that hummer al meth 


used 
A program for 
the simultaneous set of transcendental eqs ‘ 
L2 uSsIng 11), on the IBM 704 
written by the Bureau’s computes 
numerical solution was made in. the 
vested by Pickett [5 \ppropriate values OL pa- p ij 
and uw were chosen The two rea 
Wie 


ods must he 
numerical 


obtaining a solution 0) 
ri al 
COMpuULer Was 
section Tl 


Manner stg. 


values of y~ Satisty- 
were determined by successive  PProxi- 
mations These values of were then substitutes 
In eq 1 | and the values of [/2 
for each of the first three 
determined by 
rection factors, 7 
relation 


eg 7 


@d satisiving eq 1? 
Vibration wer 
approximation The cor- 
then computed from thy 


modes of 
SUCCEOESSLV¢ 


were 


Were ¢ hose hh, rang- 


Kor each 0 


values of pap I: 


1 


it) to U 


Twenty-nine 
ing in value from 1.0 
values of pup / 
from 


3958 


these eleven values Ol u wert 


chosen, 
There 
the first three 

The values of 7 


for equal intervals of 2a 


raneineg 

thus 
modes of \ ibration 

/ ind / 


Pbv a sixth-ordet 


zero to O.5O in steps ol 0.05 
were 319 computed pomnts for each o 
were then compute 


Aitken’s ‘ 


interpolation method The overall error mm thes 
computations Was less than three parts mm 10” fo 
values of pa p fC greater than about 0.0007, corre- 
sponding to d// 0.1 for the fundamental mock ig 


lower values of pa?p?//, however, machine round-off 
rror became significant In order to obtain ace 
rate values of 7 for patp? ke 0.0007, which corre- 
sponds to long thin rods, it was necessary to mak 
n first-ordet upproxim ition to the set oj Cs ry ’ 1] 
and 12 This Wis aceon plist ed iis follows 
If we let 
] ) / 
| 2 i B. rnd a? 

I: 

ther irom eas ‘ (| ) 
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} 
j 


} 


mance | and Furthermore, it is known from the elementary solu- 
v1Ve BP ai BP . BP a tion that 
= / i i 
gir pe : | 
» ¥ ) : *% ° 2 MW 
sly for 24 354 lim cot = lim ( coth = ) cot > 
thes e a as 


e fr | | 


— 
~~ 
— 
- 

aa 
— 
t 


, 5760 
het}. ' 
for odd-numbered modes, and 
ion of | Ifeqs 14), (15), and (16 are substituted in eq é 
terms ol order ¥ cancel terms of the form yp are l ’ 7 / lj { | 7 l ’ Mm 
elit . ill i il =U 
rw higher order Inspection of the lowest order terms ovinouns 2a _ — 2a - 2 
a , 
TY present shows that 
Su m ‘ ) - » ’ " 
x 2} ) 2y/ ly for even-numbered modes. We therefore make a 
; y . rm , . . . . 
A wel . first-order Taylor’s series expansion of these functions 
eal But obtaining a solution to (12 _ requires that the about the point m/2, FIVINne 
Proxl- | naconitude of y? and 4 be different We must, et = 
l { ‘a . . y / 
Cute therelore, make a second-order upproximation to ry > ; 
9 te : > un . . : Wi -fl “ 
{ (12 q (7). Also, it is obvious from eq (17) that 7? Is of cot = cot —— : : 
wer , . 2d é , } 
= on order VP, not ol order a therefore, we must keep sin- 
‘Or- " . - , . » : ») 
n tl the first three terms of (15) and (16 for this approxi- - 
le | . . “¢ 
' mation When this is done, and algebraic simplifica- : ; : 
tions are Made, we arrive at and similar approximations for the other three 
functions 
» The quantity pa-p- kK was eliminated by means of 
ae — 2 , » u a k. A : —— A 
is } 1/ » aa l ie 12(] } U, (1S 13), and the various powers of 7 which appeared 
ua were eliminated by means of a Taylor’s series expan- 
sion about the point d//=0, 7 This expansi 
where the last two terms are of higher order than the we sgt — pou « , 1 l. . Phi a 
} ! can be shown to contain no terms of order d//, so the 
rang- / first two In order to solve IS) for the two real . , 
} . first correction term in the eXpanslon Is proportional 
ich of | values of y~, assume 
to dil 
wert ’ ; . 
0 0: > ( ,P DP 19 With these substitutions and further algebraic and 
J . | trigonometric manipulation, we arrive at the follow- 
— Ss _ where ¢ and DP) are functions of u to be determined ing equation 
The substitution of (19) into (18) results in a poly- 7 ) 
utes : , “ws a m-(sin 7 sinh m) 1O+-15 | 
:  [g nomial my / ‘ If we neglect terms ol order higher ] | ( / t / )( pu a ) 
1s ; : ’ 
man p°’-, and equate the coefficients of each of sinh m—sin m AS(1 +p 


these 
he remaming powers of P separately to zero, Wwe 9 


“a rrive ul the following values for ¢ and 2D). ( vm sin = sinh = )( re - NG ), 


sin m—sinh m 


vi 
to ( ) ; . 
id-off where the upper (plus) sign refers to odd-numbered 
ace modes and the lower (minus) sign .refers to even- 
ore: dD 104-15 4° numbered modes of vibration. 
al | 6(1+- py Inserting the values of m from eq (3), (4), and (5 
Tak 


we obtain 
— “Te er a | .26225u+-0.20980u d 
4! 9 ./f 7 +-( ae : . : )( E ), (20) | 4i~1+4 SS669/ )\( / )» 


+ 


] 13.5234( l 1.3547 2u4+-0.28377 ) 4% 9° 


V, ON a ie |e 21 _ 
\ It ( (1 f )( ht ) | + 1.39469u+-0.31575u d 
(7): 


T, ~14+-26.1028{ 24) 
Equations 20) and (21) were then substituted in ‘ M l 

i —- 2), and the functions om vl <a, cot y/ <a, \ plot of T versus (d/l Is nearly a straight line. 

14 oan yi/2a, and tan y //2a were approximated in the This indieates that the first term (of order (d//)2) in 
olowing wrens “9 eng Y 7" Y iS Unagimary, 80 the Tavlor’s series eXpansion of T vives most of the 

{ following relations I old ile ritie ally contribution over the range of interest Thus. the 

s sy function 
cot. i coth — T—! 


TT 2a HH (25) 


tan — i tanh , 
Za 2a Varies ONLY very slow ly with dit 
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The method for obtaining accurate values of 7’ in 
the region (d//)* I how apparent. 
H versus (d 1); is made On a highly expanded scale 
for each of the first three modes of vibration and for 
each value of My, with the intercept being obtained 
in each case from the appropriate one of eq (22 
(23), or (24) and the other points being taken from 
the computer solution. When plotted in this manner, 
the first few points from the computer solution did 
This scattering of points 


Is 


not form a smooth curve 


was an indication of significant round-off error at 
these low values of (d//): These points were there- 
fore ignored, and the smooth curve, obtained at 


high T values of (d/l was extended to the chosen 
intercept. Finally, the values of HT for equally 
spaced intervals of (d/l were read from this smoothed 
curve, and the values of T were computed from eg 


(25 & 


4. Results and Discussion 


Tables 1, 2, and 3 give 7 as a function of d// and 
u, for the first, second, and third modes respectively 
These tables are accurate to about three parts 1n 10”, 
In order to calculate Young’s modulus from measure- 


ments of the flexural resonance frequencies ol 
cylinders, the following equations, based on eq (6 
are used: 
by | oe 
E—1.261886! ir. 
(f2 
, _ »/ = 
E—0.1660703 Fes 
TABLE l. The CO clion tacto oO the nia yf > oO 
") 
‘7 ° 
wet | sis 
{Rf {ROR Py 
» ] 4 
IRS458 ‘ g 
1084 14 
13 17S s 
{ . s] s4 . 
4:41.48 ~ — 
9 


A plot of 


Lr 


for the first, second, and third modes respectively 


\ linear interpolatior ol / versus d/l) or 7 
versus uw can be made with an accuracy of better thay 
one part 1n 10* over the range of the tables. Ty 
obtain higher accuracy Liiken’s interpolatio, 
method should be used 

In order to obtan iXlmMum accuracy irom this 


it 


resonan 


Young’s modulus determination. 
that, feasible, 


frequencies of the fundamental and first two over. 


method of 
recommended whenever 
tones be measured, and the values of be compute 
An average ot these values ol Ee should hy 
more accurate than any ol them alone The reaso) 
for this is that small irregularities in the diameter of 
a@ specimen, which ca be eliminated eo) 
pletely, will affect eac! ode differently: 1.e.. 
modes will have their frequencies shifted upward, } 
and downward The magnitude of the dis 
agreement Ie indicates tl 
overall untlormity of the dimensions of the specime! 


for each 


never 


sor 


SOTILEC 


among these values of 


The correctness of the se computations was check 
experimentally, by amet od which will be deseriby 
in a subsequent publication. Table 4 gives a com. 
parison of experimental and theoretical values of 7 
for two of the most ac urately machined specimens 
The discrepancy between theory and experiment Is 
seen to be of the order of 0.03 percent, except fo 
T The experimental value of T Wiis difficult 
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Publications of the National Bureau of Standards” 


(Including papers in outside journals) 


Infrared spectrum of hydrobromic acid, |). Kk 
ch NBS G4A, Ne ) L960 


d HBr“ has been made The two 


HBr sie apes Bt 
ive bee found to be 2558.94 ind rOs.o respectively 
Rotational and vibrational constants have been calculated 
( Tye l 1 i | ( Col " ‘ vood wree- 
: will pre VUE reported \ Lillies The centers of the two 

mol bands were used to caleulate und » and thes 
a i to be 15.58 and 0.072 respectively for 
HB ind 45.56 and 0.072 © for HBr® ; 


The spectrum of singly ionized atomic iodine (I 11), W. ¢ 
Martin and ¢ H. Corliss, /. Research NBS G4A, No. ¢ 


| [ir spectrum has been excited in electrodeless lamps 
1 photographed fr 655 to 11084 A Wavelengths and 
nated ‘ re given for almost 2400 lines \ 
si ind ¢ ens of the earhler analyses of this speectrun 
s inere d the number of known even levels from 43 to 
24 and the number of odd levels from 55 to 190. New 
gJ-facto ire Ve ce 16 levels, and the revious designa- 
s of 40 levels at hanged Improved measurements it 
I ‘ egk give 1 Col m of 7.4 © 
be subtracted f the values of all higher levels listed in 
} , da Vol LQDS ] ipproximately 
ROO « sified | = now include all of the strongest lines 
I ine S, of the ground configuration 5s?5p* is been found, 
{1 this configuration has been fitted to intermediate cou 
gy the NiIag ce dipe insite between levels of 
ground configuratior P,—'D 7282 A) and 4P,-'S 
160 A ive Dee Obst ed a 1 thei iture confirmed by 
LC 1 {Te The lin p x 1) ows hyperfine 
( ire \\ ict | ! iIpproximate igreement with i Theo- 
il calculation of the expected strueture New levels have 
: owl ries i rene extended d new ones found 
) ‘ ot e new ser “ »p ‘) d-12 G, the principal 
‘ ‘ vy it 7% 4504 a 4 has been de- 
i | sults of 1 i Ss are pared with theo- 
‘ 1) ( ~ Lr 7 CASES 


rhe third spectrum of gold (Au m1), | 


VBS GEA, Vo. 6, 481 (1960 
The spark spectrum of gold has been photographed in a helium 
sphere from 500 to 6600 A About 500 lines have beet 
g ! t hie 1 spe 1 Audi, and separated from 
si ying to « iges of ion ion, by observatiol 
po of e lines Sixt vo levels have bee 
| ls ‘ ( iris vy tro | 77, ind wi” OS CO 
Y ind 45 odd leve belong y to the 5d op ind 
6s Gp config All of the expr i levels from the 
hg 5 5 ) 6s and 5d* 6p ive been identified 
Au \\ hese levels vas possible to classify 256 lines 


Note on particle velocity in collisions aereen liquid drops 


and solids, O.G. Enge Research NB 506A, 197 (1960 
e developed for plane-way a velocity 
mauced ! Nid-against-liquid COLISIO!I \! explicit 
press | | ( ol less ( flici a | 1 ippears 
these s is deduced 


Error analysis of a standard microw ave phase shifter, G. | 
Schafer and R. W. Beatty Research NBS 64C, No 


proposed 


iched to a 


phase shifter = heer 


classes oO 


MW MM ; ; 
errors Col lered are those due to imperfect tuning of the 
¢ ie ction and those due to inaccuracies in deter- 
gy the motion of the adjustable short circuit and the 
wid of e waveguide termed tuning and dimensional 
rors respectively Cit iphs ire presented for estimating 
! itter rf coefhiel magnitudes needed to estimate 
fi tuning error from observations of amplitude changes 
! tuning procedure Graphs are presented for esti- 
mits of dimensional error for WR—90 waveguide in 
f ommended frequency range of 8.2 to 12.4 kMe/s from 
olerances of the ixial motion of the short circuit and broad 

ons of the waveguide 


Electron resistivity studies on the Athabasca Glacier, 


Alberta, 
Canada, G. V. Keller and F. C. Frischknecht, J. Research 


VBS 64D, No. 5, 489 (1960 

Phe use of electrical methods for measuring ice thickness and 
properties on t he Athabasca Glacier, Alberta ( nada has 
rent tudied by the U.S. Geological Survey Two methods 
fo easuring resistivity were tried: one, a conventional 


which current was introduced gal- 
glacier through electrodes, and the other 
method in which a wire loop laid on the 
ce was used to induce flow Results of the 


resistivits method in 
inically into the 
omagnetic 
“nic 


current gal\ 


re irements showed large variations in the resistivity of 
the ice: in a surface layer several tens of feet thick the resistiv- 
tv is between 0.3 and 1.0 megohm-meters, and under this 
ive! tin resistivity of the ice is more than 10 megohm- 


resistivity of the surface ice determined by 
vater content rather than bv molecular 


had no effect on the mutual coupling measurements 


resonnunce loss. 


n the frequenev range from 100 to 10.000 eveles per sé cond. 
As a consequence the electromagnetic data could be inter- 
preted simply terms of ice thickness and bedrock resistivity. 


Amplitude distribution for radio signals reflected by meteor 
trails, A. D. Wheelon, J. Research NBS 64D, No. 5, 449 U 


Ihe probability distribution for the envelope of the received 


At) 


ignal composed of reflections from many meteor trails is 
derived theoretically soth the effects of numerous, small 
meteors and the residual reflections from infrequent, large 
meteors are treated simultaneously For ‘the particular eX- 
imple of exponential decay of initial spikes which are them- 


then 


composite 


tributed is the 


amplitudes, 
residual 
given by 


mnverse square ot 
that the 


plitude exceeds a prescribed level r is 


probability 


l 
PCR = 


[+ a 


This functic behaves as a Ravleigh distribution for small 
implit ude irgins For the larger, less likely amplitudes 

irees with the result predicted by elementary analysis of 
olates weteor reflections Possible refinements of these 
resu ire also discussed 


Computation and measurement of the fading rate of moon- 
reflected UHF signals, 8S. J. Fricker, R. P. Ingalls, W. ¢ 
Mason, M. L. Stone, and D. W. Swift, J. Research NBS 64D, 


\y Ae) 
A id is described for pre dicting the fast fading rate of 
noon-reflected signals It is based entirely upon considera- 
0 e observer-moon positions and relative motions 
Iexperim il results which are in good agreement with the 
comput | fading rates have been obtained from a moon- 
reflectio experiment at a frequency of 412 Me/s porn 
ossible in pli itions of this method of it terpreting fading 





On the theory of wave propagation through a conce ntrically it very low frequencies, decreases 


stratified troposphere with a smooth profile, H. Bremm«: Lum ¢ i 
J. Research NBS 64D, No. 5, 46 1960 queneyv, then incrense I ") it the @g otrequeney The su 


The W.K.B. approximation for the solution of the height velocity of energy flo lep I 1] tro he lor P 
gain differential equation for a curved stratified troposp ( tudinal value excep It es OF rv sma pr 
is discussed in detail The approximatior l¢ pel ds mainly o1 wave ormal angles 

a variable wu which can be interpreted as the height de- | 
pendent contribution of the phase for a field solution obtained | Shielding of transient electromagnetic signals by a thin 

by separation of variables An expansion of with the conducting sheet, N. R. Zitron, J. Research NBS 64p, ‘ 
aid of partial integrations leads to further approximations Vo. 5, 563 (1960 

which facilitate the determination of the eigen values ind of Ihe sniesding effect of I . horizontal imperfectly CON. 

the amplitudes of the modes connected with the propagation duet sheet against the tr t field of a vertical magnet sp 
problem. The influence of the refractive-index a if as- dipole when excited by a ramp function Is investigated, Ths 3 
sumed as smooth ther appe ars to be restricted to depend- results are calculated 1) { | place transtorms of the fo 
ence on the surface values of this index and of its gradient frequency spectrum fur s for the steady-state probler 
insofar as propagation ove! the ground is concerned Further The re sponse to the I p Ttunctior is calculated and g 
all height effects of elevated antennas can be expr pea) In significance of tl re hielding igauinst surges js \{ 
terms of the distance to the corresponding radio horizor discussed 


This results in simple relations between the fields connected 

with two different refractive-index profiles, provided bot! Cylindrical antenna theory, R. H. Dunean and F. A. Hinche 

profile coincide near the earth's surface J. Research NBS 64D, N 69 (1960 fo 
\ partial survey of evlindri intenna theory pertaining t sit 


t model with a rrow gap is presented Che surve 


Polarization and depression-angle dependence of radar ter- 1 tubular 
rain return, I. Katz and L. M Spetner, / Research NBS 64D, includes discussion of the tl ( if Hallén, King and Middle- \ 


No A. iSO 1960 ton, Storm, and Zuhrt \ ptual relation between theor rl 

A study a nn rae on radar back seatter ind experiment is described he latter part of the article is ) tl 

ing from land and sea surfaces indicate i) the polarizatior concerned with a new Fourier ries solution of the Hallé st 

dependence of the normalized radar cross section, o of equatiol Mhis solution 5 a eloped suct wavy that tl 

ocean surfaces cannot be explained by the usual ‘“‘interference expansion coefficients are 1 inknowns of a svstem of line 

phenomenon’’ and (b) there is a distinct difference in the forr equatio Phe ele f coeffici matt ire gi 

of the depression-angle dependence in that o» for ‘smoot DV a ghiv converge N ‘ | results are give 

surfaces follows a negative exponential whereas oo for ‘rough for half and full wa I ‘ le vith ! I lengt! 

surfaces drop off as the sine of the depressior ingle radius ratios of 60 and 500 These results compare quit 
closely ose obt | King- Middlet eCOry Ca 


Methods of predicting the atmospheric bending of radio rays, = : ; , , : - 
B R Bean. G. D. Thayer. and B. A. Cahoon i ceil Vibration-rotation structure in absorption bands for the 9 
calibration of spectrometers from 2 to 16 microns, E, K \ 


VBS 64D, No. 5, 487 (1960 ._f . . . 
Three methods for predicting the bending of radio ravs whet Plvler, A Danti, L. R B eC, nd E 1) lidwell, NBS 1 
et refractive index profile above the surface laver is unknow1 VWono 16 (1960) 20 cer 
ha e been developed recently bv the authors These methods Suitable bands fr comn - es | deta aren t ibulated al 
are: a statistical technique for refractio it high initial ele remeasured wherev from £ to 10 Ons 

vation angles, estimation of bending from a exponential obtall in accuracy ol ag U.U0 CI throughout the regy 
model of at mosphe ric refractive index, and a modific yf ind to provide good UID point it Trequet! A ntervals 

the exponential model to account for the heavily weighted Pome GUY rotatio ' llustrated V spectro- 
effects of anomalous initial refractive index gradic ts at miall e . Nave d Dp - tables 
initial elevation angles Kach model is dependent upon the with considerable ;, rons obtan 

value of the refracti index at ground el « Be ise of nh other taborators . rpl n bands wel remeasul ' 
superrefraction, the additional knowledge of the refractive hee age spasegh nn Pies f ciyr ey, & i - era g 
index gradient next to the earth’s surface } ethod ; ' a“ . . ~ 5 =" uv 

works best in a particular range of initial elevatio ingles or ' aed. Perot I ( I terist features 
meteorological conditions. The height and angular nges of x " ind tual band = rt priel na goer De 
gee ee oe Tee ee | calihenl clude H,0, CO,, CO, HCl, HBr, NH, C)Hy > ¢ 
representative of wide climatic Variatiol It =s Tou i tf it i 4 N20 nd po Repr ted ft the Je irl _ 
the use of the best of the hree methods w ilwa - ilt . I Rese or tt! Nat - : . if Standard \. Physies 
i predictior of the total atmosphet! hend } ner 1 4 tr Vo H4.\ | ’ Februar 1OG60 

cent for initial elevation angles from zero to 10 

within 4 percent for initial elevation angl ren 17 Dynamic measurements of the magnetoelastic properties 

mr (~1 degre: of ferrites, V. KE. B VBS Tech. No + (PB151408 , 


Use of logarithmic fre que ney spacing in ionogram analysis, 


G. A. M. King, .. Research NBS 64D, No 501 (196 I 1 : ' 
The use of logarithmic frequeney spa the red p ry I 
ionograms to electron density profiles br I fund 

tages Among then s the fact t! | ! I I 

for the analvsis, one need oO tete re 1) . 
index Formulae involving only the phase ref ( f u 

ire presented; for the ordinary componet one exne ind om Cle 1 | jue ( tol 
ipproximate formula are given, while for the extraord ! ist 1u per ! tostrictior ve fficle En 
component there is at ipproxin ite for 1A " i t 1 loss Tact the fe t eter ned \pparat sol 
wide range of geomagnetic latitude | ( | f d p ) i p 32, 
cussion of quasi-longitudinal appro itions to the extra M 


dinary phase refractive inde» Magnetic drum directory and programming system for 

codesorting letter mail, I. ‘ VBS Tech, Note 4 
Guiding of whistlers in a homogeneous medium, R. |.. Smith, PB151 1960 . vil 
J. Research NBS 64D, No. 5, 505 (196 Phis report is an al 1 extension of the Rabinel r) 


The velocity of energy flow of whistle il i homogeneous Engineering Con pans ' f drum fi : 
medium is computed is a function of wave-normal angle director ind special purposs D i 1 it 1 tes ( 
Phe maximum allowable cone of rav angle ipproact yP OC) vied ddr P | ‘ , syste i 
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che 
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scheduled for lat 
Post Offices The 


examine the 


1960 installation at the Washington, D.C 
is composed of 6 major sections that 
considerations involved A 

first Then the general 
e tie lds COTN posing iddre Ssses 1s discuss d 
facilitating both the human 
computer coded address-to-bin- 
Addresses are classified into two types: 
standard When coded, their fields 
sufficiently flexible to specify a 
must furthermore adequately 
forms of misaddressed mail that 
manually \ slightly amended 
then presented and analyzed in 
specified for coded addresses in 


re port 


successively major 


short introduction is presented 


problem of coding 
This ares is al 


y process and the 


ilvzed in terms of 


imber 1 inslation 

standard ind ‘‘non 
must be nique ind vet be 
Chey 


more Common 


variety 1 
specily SOTTL 
s presently correctly sorted 
Coding Plan E-1 is 


Lhe conditions 


form ol 
terms Ol 
gene ral 


After the 


presentation Of the mail sort 
discussed 


mail 


coding area, the local 
Rabinow 
extended to 


ind 


program ot the 
detail It is 


translatior program 


proposal is 
include the 
non-standard and special address 


and it 


ther outgoing 


forms Information istruction characters neces- 
sitated bv the programs are listed in detail ind the effects 
of the file-directory source document preparation determined 
\ coding keyboard is then suggested 

The programming is then reexamined and it is determined 
that all iddresses Cal be treated is special cases of the 


standard iddre ss forms at il) INCTeCUSE d cost In memory space 


ind with a 
Finally, alter ililve 


saving in the special instruct 
methods ot 


ire ¢ xplored and 


ion symbols required 
information 
terms Of memory 


fil -directory 
Valuated in 


and 


storage 


spac 


‘ 
requirements, look-up directory ACCESS times, 


ind conceptua s 


Carrier frequency dependence of the basic transmission loss 
in tropospheric forward scatter propagation, K. A. Norton, 
NBS Tech. Note 53 (PB1614554 I 
\ further interpretation is given of certain Lincoln 


lata obtained in an experiment scaled 


Ht) S700, 
Laboratory 


using 


antennas as 


presented in a recent letter to the Proceedings of the I. R.E 
from Bolgiano Phis paper has four objectives: first, to 
irify tl significances of these data from the st indpoint of 
e engineer d veloping long-rangt tropospheric scatter 
systeln.s irther statistical inalvsis to the se 


thei irds the theory 


second, to apply af 
| significance as reg 


f radio propagatlol hrough a turbulent itmosphere and 

fourth, to describe a suitable method for the measurement of 

( I eorolog i Dar eters ¢ 1 y | neory Based 

nis a il\ 0 | Col | wmoratory | La, it is concluded 

it the carrier frequene dependence of the Dasic transmission 
ss ( ot de il ble trol nour to hou 


Determination of a general index of effort in sorting mail by 
conventional methods, S. Henig, NBS Teel Vote j 


PB161 , ] Ht) ‘/ 
leseril 1 OV nro I reference Oo iVailable 
statist il 0 of letters origi Ing Within thie 
f S st i bv s irge post offices will be involved it 
I ! in five iL to be compte el ted to a carrier 
f ( of boxes \ met 1 presented for the 
‘ . ( 1 du ! lings which are 
lefined of measurement fe dex of sorting 
Tor It found that the ippel bound erage number of 
lings pr ‘ ic ree large post offices varies from 3.16 
2.98 and tha bounds are applicable to 96,2 to 92.4 
f il originating first class letters t is inticipated 
it the values of such indexes for all large post offices will bi 
nsistently close to the range alread letermined 


of N., O., and NO molecules trapped in 
Broida and M,. Peyron, J. Che Pp} 


Emission spectra 
solid matrices, I 


32, Vo _ 1068 rT, 

Moleeular systems are observed in emission it solid products 
Irom a on iiscl 0 ipped i lig 1 hel temperature 
Pre Ol 7 ilar issig ne! ~ ‘ Der checked 
vith t elp of isotopic substitutions of oxvgen and nitroge 
The Herzberg svstem (A bands) of ox 18y X % 

Ss analyzed and the molecular constants are derived for a 
molecule trapped i nitrogen matriy Another svsten V/ 
Dands) is attributed to the NO molecule (‘II V?1l 


expressed in terms of the 


The components of power appearing in the harmonic analysis 
of a stationary process, M. M. Siddiqui, presented at Symp. on 
Statistica Vethods in Radio Wave P 
Los Angeles, June 18-20, 
Statistica Vethods in Radio Wave P 
Press, Ne York, N.Y., 1960 
Suppose that from a continuous record of a stationary process 
over a time interval 7’ seconds, 2n l equally spaced readings 
and the usual harmonic performed 
Assuming the true frequencies present in the process to be 
Li/3T). WaT), Wz. Bz. evcles per second, 
that the power ascribed to the frequency T e¢/s 
consists of three components: (1) the true power in the 
frequeney j/ 7 ¢/s; (2) the powers in the frequencies (2kn + j)/ 


ropagation, Un ersily of 
1958 


Reprinted 


opagation, Pe 


(a fornia, irom 


gamon 


ure iken analysis is 


it is show! 


3, / 3) a part of the powers in the frequencies 
| Tr’) e/s,? 2,3 
Introduction—Can you measure it? L. V. Judson, Am, Soc, 


lool & Mfgr. Engrs. Tech Paper 239, Book 1, 60, 1 (1960 


\ papel which is the introductory paper In a Metrology 
Symposium It indicates the need for precise and exact 
measurements and it treats two subsidiary questions: (1) 


and 2 


\W hv measure ? 


ment? 


What has been the historv of measure- 


Correction for systematic wavelength shift in atomic beam 
devices, R. L. Barger and K. G. Kessler, J/. Opt. Soc. Am, 50, 
Vo. 4, 352 (1960 

Light emitted or absorbed by atomic beam is 
shifted in frequency relative to the frequency for a stationary 
itom if the light ray is not normal to the trajectory of the 
m The magnitude of this shift is calculated and several 
which correct for this effect are described For 
recording of interference fringes, complete 
achieved In the case of photographic 
that the residual error is negligible 


atoms in an 


} 


aevices 
photor lectric 
can be 


shown 


COTM pM nsation 


recording, it is 


The error in prediction of F2 maximum usable frequencies 
by world maps based on sunspot number, I. L. Crow and 
| ) H Z icharisen, Seclion fe Radio Propagat on Pi enomenolo- 
i, book, Statistical methods in radio wave 
Pergamon P Vew Yo - \ Pas 1960 
The feasibility of preparing semi-permanent world maps for 
predicting /2 maximum usable frequencies based on sunspot 
number is studied. The components of variance of the predic- 


p opagal on, Pp. 7i8 


tion are estimated; they arise from the following four major 
sources of error: determining regressions of /'2 characteristics 
on s Inspot number with limited vears of data, the seatter of 
Do 1 about these re LTEeCSSIONS, drawing world contoul maps 
with limited numbers of stations, predicting the future sunspot 
number All four sources are found to contribute about 
equally to a total standard deviation of monthly median 

iximum usable frequency predicted five months in advance 
of about 1 and 3 Me for 0 and 4000 km respectively Hence 
the ] posed maps are concluded to be feasible 


Bremsstrahlung linear polarization, J. \ 


Placious, /1 Nuovo Cimento 15, Series \, 571 (1960 
Tl ork presents a general quantitative deserption of 
1e bremsstrahlung linear polarization on the basis of original 


theoretical caleulations 
he 


experimental lata and available 


give the deper lence of the polari ation on a) t 


nitial electron kinetic energy J) in a range from 10 to 
10°? MeV, 6) the photon energy in a range from 0.1 7» to 79, 

the photon emission angle in a range from 0 to 180 degrees, 
ma the atomic number of the target in a range from 4 to 
79 The experimental data were obtained for a range of 
electron energies from 0.05 to 1.0 MeV with beryllium, 
iluminum, and gold targets Theoretical estimates of the 
polar itlol were calculated from the Sommerfeld-Kirk- 


patrick-W iedmant results tor the non-relativistic 
region, the Olsen-Maximon results for the extreme relativistic 
on, and the Gluckstern-Hull (Born approximation 
ntermediate regiol Final results are 
peak polarization and the 
a function of the electron and photon 
polarization are 
basis of the combined « xperime ntal and theroretik 


energy 


energy 
corre- 
ingle as 
given on 


il data. 


best estimates of the 





Geomagnetic disturbances and velocity of slow-drift solar The wavefront reversing interferometer, J. B. Saunders 


radio bursts, M. G. Wood and C. S. Warwick, Nature 184, ( ” on O7 Sponsored he Int é 
No. 4697, 1471 (1959 ( of Op } / V/ . 1958 
Frequency drift rates were measured on record ) e Ke P (/p Vie j / \, 
Davis dynamie spectrum analvzer for 48 Type II (slow-drift } V.) of 
solar radio noise bursts The change vith trequeneyv ot the Phe Wwavelrol! revel ‘ ( 1 lO tes 
frequency drift rate is found to differ for two groups of burst the shape of wavetront If ( 
those followed by geomagnetic d rbance and those followed e@ sfhiape l Opti ( ! 1e1 
by no disturbanee This change in drift rate is interpreted as he evaluation of the ‘ il! , ition of 
acceleration of the radio soures is it moves outward throug the element It serve to easure the shape ‘ flecting 
the solar corona Bursts followed by geomagnetic disturb rraces ag iberratlo ro 1 ons of 
ance determine i rate that correspo is to enute me ( Opt | pa ns il i beat mogenell I ( medi Mi 
acceleration of the souree throug viiict . - I \ ghee ethod « 
inaly gy the data - \ Str ble f | 
It t l ertere 
The outlook for machine translation, F. L. Alt, Proc. VW ng ( be « ed bv 
Joint Compu Cont lo / D N } vudi« ‘ ‘ | \ 
Calif., (May 1960 17 
The current status of the maior research proiects concerned Spectra emitted from rare gas-oxygen solids during electron St 
with the use of computing machines for the translation of bombardment, L. J. Sc! H. P.B la. J. CI P} 
natural languages (e.g., Russian to Englis! s surveved $2, \ 3 Lost 
None of these projects has as vet reached the stage ot Tull grit ¢ s10n Trom rat i I LiTiigis mour 
automatic translation of satisfactory qua The difficulties | Of oxygen has been « electron bombardme 
encountered fall into two man classes <emanty nd technique Phe Her | I U 1 e ore \ 
tactic Several methods of ipproach have be sed elieved ! the ») } I ‘ OXV ge 
overcome them; these are surveyed 1 the approac f the te fo ( obse 1 sne . 
National Bureau of Standards is described ! more deta I i! I l X, Uh ~ i N ‘ ot ©) 
It appears likely that within a few vears automatic translation Heel lentined The eff I é Is 
of fair quality may be possible it a cost substantially below prationa tru re ol H perg band CIs ed 
that of human translatior The m yor cost element nitin , 1) 
development of machine programs, cost of macl “¥ Accurate microwave wavemeters with convenient calibration 
operation—are analyzed in their dependence on th athe tables, H. IX Busse iA. 3. 8 I? S I) 31 . 
of translation and on the tvpe of equy ent used If 4 , 
Ler rs 
Statistical evaluation of interlaboratory cement tests, J. I s. High Q's wer efined tructi 
Crandall and R. L. Blaine, Am. Soc. Testing Materials Pro technique \t 9000 M f 0.02 Me was obtained 
59, 1129 1QG59 | ( it) Llite CCU ! | I rains be 
A cooperative series of phvsical and chemical tests was mad the meta ( \ . 
, : . lt)? entre is for / f } } rhe } mt 
bv 103 laboratories on twelve samples of cement during the p 
period of one year Results of tests bv all laboratories o1 ar ; ‘ re met ist e fitting 
each property were plotted | sentter | oran iccordnu : : ‘ ' v ribed 


the Youden method to enable the participating laborator » 
The Stokes flow problem for a class of axially symmetric 


to quickly evaluate their results Che 1 J . 

| > | " 17 - ) 
St itistiecally Bot} Statistical and raphy met wis vere bodies, | | W H / ' , l 
emploved to indicate the precision of the test methods and the fe ' Jo 
extent of laboratorv bias The lack of precision th some The Stokes flow ] ' agees 19 
of the tests indicated a need for improved eat a fF goo pr wagtiic: 
Tr \ f b, f +} ] , 

‘ Oo of each of é wus labo ‘ f 
by a rating svs mm A large 1 ry | ( ‘ ot f D 

| ’ i?’ ‘ i | i ‘ ‘ ) 
showed poor agreement with the larger roup Niar of the F ; , 
' i ? i} 
liscrepant laboratories did not Mmpro iu! the one ] { rer os . ] 
period the test program Was active r ; ae } j | ps 
Discrepant results optaine | 1) a Te I f f \ \\ ; \ 
ippreciably trie Stal lard 101 tions trong l ol 
The various tests thre ’ i ‘ el ! ‘ = - 7 
: 7 ) 
indicated 
r ~ ) j a 


Effect of notch geometry on tensile properties of annealed nit 
titanium at 100°, 25 78° and 196° C, G. W. Ge ind hod J 
N. L. Carwile, Am. Soc. Testing Materials Proc. 69, 985 (1959 Exnlicit ca a ee so ' oy _ 


a circumferential notch on the 1 e |} , ( t \s ex- ; 
annealed commercially pure titanium Wa 1! est | 1 Ihe I ! pr 

notch geometry was varied by changes in its dep nd ot 

radius; the notch depth, in terms of the cross-s 

removed in machi o the ote} ' if a QJ p 


and the root radius ranged from 0.005 to 2.0 ir Prue re Analysis of two-factor classifications with respect to life tests, 


true strall data were optaine i to the mitiatior it | ture \I if el hao ( j \ 

of specimens at.100 aD”, 78 ana 106 ( nd repre j pe j US >) } / ‘ ( 1960 ¢ 
sentative true stress-true stra curves are pre nted ‘ Recent i ‘ of t! 43 
relationships of selected strength and ductility indices to the | expo listril pp \ 
noteh geometry and to the iccompal ! tria i na | f » (yy 

stress concentration factors ire presented graph nad Sob I} ! i ! es Hy 
discussed in some detail in the text procedur ‘ | e ¢ 
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onential dist! Dil 0 This papel outlines procedures to tx 
ised for ilvzing the results of a two-way classification wit! 
spe tine scale factor of the exponential distribution 
The test procedures ire based on likelihood ratio tests and 
ive the al ilogs in the isual analysis-of-variance tests for 
iin effects and interac ns \pproximations are given for 
sts based on small samples These ime procedures carry 
yer directly to a two-wa inalysis of variance, if the quar 

ties 1 iy iwnalyzed ire Variance estimate following ch 


) 
rare distributions 


Microwave spectrum of trans-crotononitrile, V. W. Laurie, 
Chem, Phys. 32, No. 5, 1588 (1960 


The microwave spectrum of ethyl evanide has been studied 
the region from 17-36 kMe Both parallel and perpen- 
licular transitions have been assigned Rotat onal constants 
\le for the ground vibratior il state ire a 27 6603 30, ‘y 
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